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EVIDENCE OF THE NATURE OF ROTATING LIQUID HELIUM* 


J.R. Pellam 
California Institute of Technology, Pasadena, California 
(Received August 18, 1960) 


Earlier experiments by Hall’ and also by Walm- 
sley and Lane? on the angular momentum content 
of rotating liquid helium II have involved meas- 
urements of torque required to set a vessel of 
liquid into equilibrium rotation. Integrating the 
resultant torque-time measurements over the 
angular acceleration (or deceleration) period 
yielded thereby a measure of the angular mo- 
mentum acquired by the helium. Since further- 
more an integration over the volume of the vessel 
was involved, such methods produced composite, 
over-all data. 

The present method® permits investigating both 
the time and position dependence of liquid helium 
II circulation within a rotating vessel by suspend- 
ing a movable Rayleigh disk device within such a 
system. Besides thus providing both instantaneous 
and local observations, the sensitivity of such a 
probe to individual fluid component motions per- 
mits in fact a still further resolution; the special 
case of zero response, for example, uniquely 
specifies each fluid component as stationary. 

The liquid helium II sample investigated resides 
within a cylindrical shaped glass vessel (just 
under 5 cm inner diameter) rotated axially about 
the vertical; a small ($ mm) aperture at the lower 
(tip) extremity connects the region with the outer 
helium bath for filling purposes, etc. A “sta- 
tionary support” is provided inside the container 
for the disk suspension system which, although 
not rotating, is equipped (rack and pinion) for 
position adjustment along the vessel radius. The 
disk probe (3-mm diameter) is suspended by a 
quartz fiber (torsion constant ~ 10-5 dyne-cm/deg) 


to present an angle of 45° to the rotational flow 
direction. Resultant angular deflections from 
equilibrium orientation (tending to turn the disk 
across the flow) are observed as displacements 

of a narrow light beam refected from the (silvered) 
disk surface and focused upon an optical scale. 

The present observations are for a vessel rota- 
tion rate of 2 rpm. 

The results of this investigation include the 
following principal findings: 

(a) A delay interval of several minutes elapses 
between setting the container suddenly into rota- 
tion from rest and initial response of the disk to 
the liquid helium II system. In contrast, the disk 
responds immediately to stopping rotation of the 
vessel. Thus liquid helium II appears to acquire 
a “rigidity” in rotation which is absent when sta- 
tionary. 

(b) A linear relationship between such delay 
interval and radial position of the disk indicates 
that the circulatory motion propagates as a “ro- 
tation front” slowly inward from the moving wall 
of the container with a definite velocity. 

(c) Steady-state angular deflections, reached 
after roughly one-half hour buildup time from __ 
initial response, obey a parabolic dependence 
upon radial distance from the axis. 

(d) Finally, steady-state deflections depend 
markedly upon temperature. Upon cooling through 
the A point the apparent rotational flow drops 
abruptly from the helium I value to zero, indica- 
ting no circulation at all. With reduction of tem- 
perature further into the helium II region circu- 
lation reappears, however, and builds up mono- 
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FIG. 1. Time delay versus disk location. The time 
interval (min) preceding initial response is plotted 
versus distance (cm) from rotation axis to midpoint of 
disk. Note location of intercept. 


tonically to the helium I value toward the lower 
temperature limit (7 — 1°K). 

The “rotation front” propagation phenomenon 
(b) is illustrated in Fig. 1, where response delay 
times (min) are plotted versus radial distance 
(cm) to the disk midpoint from the rotation axis. 
For a temperature of 2.0°K (and vessel rotation 
rate of 2 rpm), the observed slope indicates an 
apparent “rotation front” velocity radially inward 
of about 2 mm/sec. Furthermore, observation 
that the intercept (zero time lapse) falls inside 
the liquid helium region indicates the liquid “cir- 
culation front” first disturbs the extreme edge of 
the disk. 

The unexpected temperature dependence (d) is 
illustrated in Fig. 2, where fluid circulation is 
plotted versus temperature (°K) for a fixed disk 
location (1.84 cm from axis of rotation). Observed 
results are expressed most directly as an effec- 
tive density ratio (pef¢/,), representing the frac- 
tion of fluid evidently carried into rotation with 
the vessel; the ratio measures the observed torque 
normalized to the helium I value (at the \ point). 
Clearly a marked dependence of temperature 
exists, characterized by a sharp drop evidently 
to zero circulation immediately below the A point. 
A general (approximate) resemblance may be 
noted between this curve for effective density ratic 
and the well-known temperature dependence of 
superfluid concentration (p,/p). Such behavior 
might be interpreted (but possibly oversimplified) 
to indicate that under these circumstances super- 
fluid rotates whereas normal fluid does not. 
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TEMPERATURE (°K)—— 


FIG. 2. Temperature dependence of fluid rotation. 
The effective density ratio (p.¢¢/p,), plotted versus 
temperature (°K) for a fixed disk location (1. 84 cm 
from axis), represents the fraction of fluid carried 
into rotation with vessel. Note apparent rejection of 
fluid rotation upon cooling through A point. 


Logical explanations for findings (a), (b), and 
(c) follow directly in principle from an applica- 
tion* of the Feynman concept of rotating liquid 
helium II, based on migrating quantized vortex 
distributions. Observed temperature dependence 
of finding (d) appears unaccountable, however, 
on the basis of either theory or (previous) experi- 
ment. According to Feynman’s theory’ and Hall’s 
measurements,‘ for example, indicating rotation 
of both fluid components with a containing vessel, 
no temperature effect might be expected. Accord- 
ing to Walmsley and Lane’s measurements, ” on 
the other hand, indicating rotation of normal 
fluid alone, an increase of circulation with tem- 
perature might be anticipated. The actual ob- 
served behavior remains clearly oblique to such 
expectations. The possibility naturally exists at 
this juncture that some obscure phenomenon re- 
lated to the method of measurement is responsible. 

A paper now in preparation to describe fully 
the experiment and results will be submitted 
shortly. 





*Supported in part by the Alfred P. Sloan Foundation, 
Inc., and the National Science Foundation. 

‘H, F. Hall, Trans. Roy. Soc. (London) 250, 359 
(1957). 

*R, H. Walmsley and C. T. Lane, Phys. Rev. 112, 
1041 (1958). 
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3The technique as such was described at the Hawaii 
American Physical Society meeting, prior to present 
findings [J. R. Pellam, Bull. Am. Phys. Soc. 4, 370 
(1959)]. 


‘R, P. Feynman (private communication). 

*R, P, Feynman, Progress in Low-Temperature 
Physics, edited by J. C. Gorter (North-Holland Pub- 
lishing Company, Amsterdam, 1955), Chap. II, Vol. 1. 








FIRST-ORDER TERRESTRIAL ETHER DRIFT EXPERIMENT USING THE MOSSBAUER RADIATION 


Martin Ruderfer 
Dimensions, Incorporated, Brooklyn, New York 
(Received June 3, 1960) 


Despite the general acceptance of the special 
theory of relativity, the Fitzgerald-Lorentz con- 
traction theory and its ability to account for 
previous experiments have persisted. In particu- 
lar, the extension of the Fitzgerald-Lorentz con- 
traction theory by Ives! appears to bring the 
contraction theory on a par with both the special 
and general theories of relativity. It has been 
pointed out? that this impasse may be resolved by 
astronomical observations and satellite experi- 
ments designed to detect a variation in the one- 
way velocity of light with direction in space. The 
existence of such an observed variation, which 
may be considered to be synonymous with the 
existence of an ether, is allowed by the Fitzgerald- 
Lorentz contraction theory but not by the special 
theory. This note proposes a terrestrial experi- 
ment which has considerable sensitivity when the 
Mossbauer radiation is used and which permits 
a first-order test for a variation in the velocity 
of light with direction in space. 

An electromagnetic radiator having a highly 
stable frequency, a frequency-selective absorber 
or filter, and a counter or other suitable detector 
are placed on a turntable, as shown in Fig. 1. 
Electrical connections to the counter or detector 
may be made through slip rings. If there is an 
ether and the velocity of the plane of the turn- 
table through the ether is », then the time taken 
for the radiation to travel the distance s between 
radiator and absorber is 

s ‘ 

"= ory 088 ’ (i) 
where c is the velocity of the radiation in the 
medium between radiator and absorber. When 
the table is rotated there is a variation in this 
time if v#0. As the time varies there is a 
change, -fdr, in the phase of the radiation as 
received at the absorber, where f is the fre- 
quency of the radiation. The time rate of change 
of the phase corresponds to a change in frequency. 


COUNTER 


TURNTABLE- 








FIG. 1. Arrangement for the detection of a one-way 
variation in the observed velocity of light. 


The relative change in frequency at the absorber, 
when v <c, is then, to first order, 


Af _ dt_wsv sind v ds P 1 ds 

SF A cB 
where w is d@/dt, the angular velocity of the 
turntable. The last term is the Doppler shift. It 
and the second term may be minimized by using 
a rigid construction to insure that ds/dt is sub- 
stantially zero. 

If the linewidth of 1 part in 10"* achievable with 
the Méssbauer radiation is used for detecting a 
minimum change in frequency, s is made 1/10* 
km and is constant, and w is 607 rad/sec, cor- 
responding to 1800 rpm, then an ether drift of 
about 5 km/sec is measurable. With a “least 
count” sensitivity of 3 parts in 10'* estimated for 
the Méssbauer radiation‘ and the above values 
for s and w, an ether drift of about 0.0015 km/sec 
is measurable. 

The axis of the table may be located anywhere 


(2) 
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along the measurement axis to permit a rotation- 
ally balanced structure. The table may be placed 
on gimbals to allow adjustment of the angle be- 
tween the plane of the table and v. To remove an 
ambiguity in the direction of v, which results 


when the absorber and radiator are nominally 
resonant at the same frequency, their resonant 


frequencies should be separated. It is expedient 
to make the separation such that maximum sensi- 
tivity is obtained, i.e., when the resonant fre- 
quency of radiator or absorber is along the linear 
portion of the resonance curve of the other. The 
separation may be accomplished, for example, 
with a Doppler shift, by mounting the radiator 

or absorber on a loudspeaker movement or simi- 
lar transducer and providing a sawtooth motion 
of proper amplitude and phase and of the same 
period as that of the table rotation. A difference 
in temperature between radiator and absorber® 
may also be used to separate the resonant fre- 
quencies. 

Since this experiment is a first-order test, 
second-order effects are neglected. These in- 
clude the transverse Doppler shift, Lorentzian 
variations in mass, length, and time, the bend- 
ing and frequency shift of light in a gravitational 
field, and any anisotropy of space of gravitational 
origin. Only if a negative first-order effect is 
obtained will it be reasonable to use the proposed 
experimental arrangement to search for second- 
order effects. 

Parenthetically, the resemblance between the 
proposed experiment and the terrestrial meas- 
urement® of the red shift in an accelerated sys- 
tem is superficial. In the latter, the counter is 
fixed near the periphery of the table and measures 
the difference in frequency between radiator and 
absorber for only one orientation of the table in 
space. If there is an ether, a frequency shift 
would result from the daily rotation of this partic- 
ular orientation, but it would be very small, and 
probably immeasurable, due to the low angular 
velocity of the earth. If the radiation between 
radiator and absorber were not perfectly parallel, 
it would, in the presence of an ether, increase 
the error in the red shift measurement, but the 
effect would be small and difficult to interpret. 
To detect sensibly a variation in the frequency 
difference between radiator and absorber result- 








ing from an anisotropy in the observed velocity 
of light it is expedient to monitor the difference 
frequency as the orientation of the radiator -ab- 
sorber axis changes. To accomplish this in lieu 
of rotating the counter it is feasible to add sta- 
tionary counters around the periphery of the table 
and to compare their outputs simultaneously. 
This experiment, because it is a one-way test, 
is uniquely capable of differentiating between the 
special and contraction theories. In the contrac- 
tion theory, the observed velocity of light, to the 
first order, is the denominator of Eq. (1), i-e., 
c-vcos@. The special theory is correct and 
there is no Lorentz (optical) ether only if v =0. 
If » #0, the observed velocity of light varies with 
angle @ in contradiction with the second postulate 
of the special theory, only the contraction theory 
can be correct, and a Lorentz ether (which should 
not be confused with a modern ether attributable 
to a gravitational anisotropy) then physically 
exists. All previous experiments have been 
measurements of a two-way light signal or the 
equivalent. The measured two-way velocity in 
the contraction theory, to first order, is the aver- 
age of c+vicos@|, or just c, as predicted by the 
special theory. The cancellation of the first- 
order terms has always produced an identical 
result (to the second order) for the contraction 
and special theories. Thus, the two-way experi- 
ments, as has been intermittently noted in the 
past, are not capable of measuring v or of de- 
tecting a Lorentz ether. A one-way test, which 
has never been knowingly performed on earth | 
heretofore, may thus be termed a crucial experi- 
ment for it provides a final resolution of the issue 
between the relativity and contraction theories. 





'H. E. Ives, J. Opt. Soc. Am. 29, 183 (1939); 38, | 
413 (1948); Phil. Mag. 36, 392 (1945). 

2M. Ruderfer, Proc. Inst. Radio Engrs. (to be pub- 
lished). 

3R. V. Pound and G. A. Rebka, Jr., Phys. Rev. 
Letters 3, 554 (1959). } 

‘R. V. Pound, quoted by W. E. Kock, Science 131, 
1588 (1960). 

5R. V. Pound and G. A. Rebka, Jr., Phys. Rev. 
Letters 4, 274 (1960). 

°H, J. Hay, J. P. Schiffer, T. E. Cranshaw, and | 
P. A. Egelstaff, Phys. Rev. Letters 4, 165 (1960). 





ible 


st, 
he 
C= 
the 


ith 
ite 


ry 
yuld 
le 


ver- , 
he 


sue 


)- 


VoLUME 5, NUMBER 5 


PHYSICAL REVIEW LETTERS 


SEPTEMBER 1, 1960 





COLLECTIVE OSCILLATIONS AND GIANT DENSITY FLUCTUATIONS* 


R. Brout 
Los Alamos Scientific Laboratory, University of California, Los Alamos, New Mexico, 
and Department of Physics and Laboratory of Atomic and Solid-State Physics, 
Cornell University, Ithaca, New York 
(Received July 29, 1960) 


It is the purpose of this note to interpret an in- 


stability in compressional modes found by Glass- 


gold, Heckrotte, and Watson (GHW).* We show 
that their instability is related to the Hartree- 
Fock (HF) treatment of many-body systems by 
Overhauser,” who claims that HF theory always 
predicts a ground state in a Fermi gas with a 
stationary density fluctuation of a given wave- 
length present. Kohn and Nettel (KN)* give a de- 
monstration that for three dimensions this can- 
not in general be the case for weak interactions. 
The present work is in support of KN. However, 
in addition to their negative result, we find that 
in a linearized theory, a sufficiently strong at- 
tractive interaction gives rise to stable HF solu- 
tions with a built-in density ripple, as well as a 
GHW blow-up. For either repulsive or weak at- 
tractive interactions there is neither instability 
in GHW nor a self-consistent density ripple. 

This is in complete analogy to the elegant anal- 
ysis of Wolff* on the spin density case. In nei- 
ther of these cases is one involved in a BCS type 
transition which takes place for infinitesimal 
attractive interactions.® Further the present 
method seems inapplicable to study the well- 
known lattice type transition of the low-density 
electron gas. A possible application, discussed 
at the end of this note, is the theory of freezing 
of He®. 

Our starting point is Overhauser’s SCF anal- 
ysis.° Since in this work we discuss density 
fluctuations we shall, as do KN, confine our- 
selves to the Hartree field alone. This is per- 
missible if spin is unimportant and brings out 
the qualitative aspects of the situation. The 
Schrodinger equation for a state k with energy 
Ey and wave function oy (normalized in a box 2) 
in Hartree theory is 


[ -(h?/2m)¥? + V(@)] og- Ege, (1) 
V@=D fog FE - Fogo. 2) 
K 


Here v(r) is a two-body potential assumed to 
have a range. With Overhauser we search for 


density fluctuations of wave number q; i.e., V(r) 
has only a single Fourier component. In pertur- 


bation theory, we have 
oe= ana ier, [et/welR) Jet 
kK q 
Le /o_ @)e*O* ’ (3) 
wa(k) = [(K+ q)? - 7] /2m, (4) 
f= fet DFG, Rip, 


= fe*!8* var. (5) 


Substituting the solution (3) with the definitions 
(4) and (5) into Eq. (2) and carrying linear terms 
in g*(k) only (which is consistent with our gen- 
eral linear development), we have 


g*=[v@)/2) = [e/w_=(k’) +8*/w-lk’)], (6) 


Multiplying Eq. (6) by 1/w@(k), combining with 
the equation for g~(k), and summing on k gives 
the necessary condition that a nontrivial solu- 

tion exists: 


1= WO) Z L/woz +1/o @) (D) 


The sum on k is for |ki<kp; |k+q|>kp; kp 


= Fermi momentum; v(q) = fo(met4Fasy, With 
this remark Eq. (7) is transformed to 


1=[v@/2 De {n+ q) - n(k)] /w=(®)} (8) 


which is the same as KN Eq. (12). 

Before we discuss Eq. (8), we first summarize 
the theory of collective oscillations in random 
phase approximation (RPA). A convenient ap- 
proach is the dielectric constant method.” In 
response to a driving force of time dependence 
eiwte-€ltl, the “dielectric constant,” defined as 
the ratio of induced density fluctuation to driving 
force, in RPA is 


n(k+q) - n(k) 


€q(w) = 1+ 4nax(w)=1 eaaly wah) -orie (9) 


The root of egw) gives the frequency w(q) for 
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which waves propagate freely in the system. In 
RPA, these, generally, will not be damped. 

For q<<kp, the function 4narg(w) takes on the 
simple well-known form 


4na-(w) = -v@ate,]} +(&/2) In i 


ni 
1+ - Fel-ove)], 


(10) 


where vp=kp/m. The term [(ni/2)é] arises for 
—<1 in the limit e-—0, hence the step function 

[1 - 6(é)]; gle) =[density of states at the Fermi 
surface|/Q; §=(w/vpq). In Fig. 1, the function 
Re[410G(E)/v(q)e(€r)] is plotted for real &. In 
the figure three situations are indicated. For 
repulsive forces, there is always a solution with 
Reé>1; Imé=0. These are plasmons for the 
case of Coulomb forces and “zero sound” for the 
case of forces with range. For attractive inter- 
actions with g(e,)|v(q)|<1, no root exists and no 
free oscillation occurs. For attractive forces 
with g(er) |v(q)|>1, if a solution exists it would 
have Reé<1 in which case we have to include the 
imaginary term in Eq. (10). In this case the so- 
lution in both the above examples is easily seen 
to be Ret =0; Imé=[1+1/gv]. In particular the 
sign of Imé changes as gv passes through -1 
with the stable side being for weak interactions 
(givi<1). It is thus found that a critical attrac- 
tive force exists (|v(q)|=[g({e,)]~*) in order for 
instability to occur. This is identical to the case 
of exchange interactions where a critical ferro- 
magnetic exchange force was necessary for a 
ferromagnetic instability.‘ 

Comparison with Eq. (8) then shows that it can- 
not be a consistent equation unless g{€ ,)v(q) = -1. 
For v(q) repulsive or g(€,)|v(q)|<1, the only so- 
lution of Eq. (6) is g*(k)=0. For g(ep)v(q)<-1, 
a Hartree solution with a density ripple occurs. 
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FIG. 1. Plot of polarizability in dimensionless units 
vs frequency in dimensionless units. 
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It is manifested by a blow-up in zero sound about 
the spherical Fermi sea and according to Over- 
hauser® has a lower energy than the latter. 

A possible application of the above is to note 
that when He’ is put under pressure, there is a 
certain critical pressure (~30 atm) for which it 
freezes. Now He’ is rather an open liquid with 
[interparticle distance]>> [position in minimum 
of v(r)]. A careful calculation (using “t” matri- 
ces for v matrices) should then show how v(q) 
grows with pressure. Presumably this will oc- 
cur fastest for q~[interparticle distance]~*. At 
the point of instability, the gas-like quasi-par- 
ticle picture will break down and freezing will 
set in. 

One last comment on the mathematics of the 
instability is in order. There are cases when 
collective modes have complex frequencies with 
imaginary parts which damp the oscillation even 
in RPA. A case in point is phonons in metals.® 
Here for small g the phonon is given by the roots 
of 


€g() =1- wy /w - draw), (11) 


where w,,”=4nne?/M; M=ion mass; n= ion den- 
sity; a5 (w) =electron polarizability. One finds 
a root at Rew=cq, where c =(m/3M)”* up << vp. 
Hence this root has §<<1. Nevertheless it is 
damped rather than divergent. This is because 
the sign of the imaginary part is determined by 
the derivative of €q(~) at the root® which in the 
present case is determined by the w)y?/w* term 
in (11). In the previous case it was determined 
by a=(w). 

Finally, it should be pointed out that Thouless” 
and Pomeranchuk"™ have given general criteria 
for the stability of HF solutions which are no 
doubt related to the topics discussed here. 





*Work performed under the auspices of the U. S. 
Atomic Energy Commission. 

‘A, E. Glassgold, W. Heckrotte, and K. M. Watson, 
Ann. Phys. 6, 1 (1959). 

2A. W. Overhauser, Phys. Rev. Letters 4, 415 
(1960). 

53w. Kohn and S. J. Nettel, Phys. Rev. Letters 5, 
8 (1960). 

‘Pp. A. Wolff (to be published). 

5L. Van Hove, Physica 25, 849 (1959). 

SA. W. Overhauser, Phys. Rev. Letters 3, 414 
(1959). 

"P. Nozieres and D. Pines, Nuovo cimento 9, 470 
(1958). 

8T. Staver, Princeton thesis, 1952 (unpublished). 
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ABSORPTION SPECTRA AND ZEEMAN EFFECT 
OF COPPER AND ZINC IMPURITIES IN GERMANIUM* 


P. Fisher and H. Y. Fan 
Department of Physics, Purdue University, Lafayette, Indiana 
(Received July 29, 1960) 


The absorption spectra of group V and group III 
impurities in germanium have been reported 
previously.» For these impurities of small ion- 
ization energy, theoretical treatment based on 
the effective-mass formalism should be a good 
approximation, and the experimental results, 
especially in the case of donors, :bear this out. 
Experimental investigation has been now extended 
to two acceptor impurities, Cu and Zn, of higher 
ionization energies. The ionization absorption 
of Cu and Zn in germanium has been seen be- 
fore*»* but no observation of the excitation lines 
has been reported. Each of these impurities can 
bind more than one hole. Studies have been made 
on the neutral and the singly charged states. It 
may be expected that the effective-mass formalism 
would not apply for the ground state but may still 
be valid for the excited states. The absorption 
spectra provide the means to check this assump- 
tion and to reveal effects due to differences in 
the ground states. 

Measurements were made at liquid helium 
temperature. For the observation of the Cu™ 
and Zn” absorptions, samples suitably compen- 
sated by Sb impurity were used. The solid curve 
in Fig. 1 shows the absorption spectrum of Cu. 
The observed lines are labeled according to the 
designation used for the spectra of group III 
impurities.? As seen in Table I, the energy 
spacings between the D line and the various other 


lines are very close to the corresponding spacings 
in the group III spectra. This applies also to the 
spectrum observed for Zn, the data for which 

are included in Table I. In the case of Zn, addi- 
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FIG. 1. Absorption spectrum and transverse Zee- 
man effect of Cu impurity in germanium. Copper con- 
centration is 1.2x10'® cm=*. Magnetic field is parallel 
to [111] in the (110) plane. 


Table I. Ionization values and energy spacings between the excitation lines for various impurities in german- 


ium, in units of 107° ev. 








Group II 
(average values) Cu Zn Zn- 
B-D 1.388 1.398 + 0.020 1.382 + 0.051 eee 
C-D 0.722 0.745 + 0.016 0.729 + 0.033 3.33 + 0.06 
D-G 1.73 1.72 +0.03 1.84 20.03 aoe 
Ionization Optical 42.8 32.6 85.8 
energies Potemens 40 30 90 
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FIG. 2. Absorption spectrum and transverse Zee- 
man effect of Zn” impurity in germanium. The sample 
is compensated with estimated 3.810'5 cm™ Zn and 
2.2x10'5 cm™ Sb. Magnetic field is parallel to [111] 
in the (110) plane. 


tional structure was observed on the long-wave- 
length side of the G line, the explanation for 
which is not yet clear. The solid curve in Fig. 2 
shows the spectrum for Zn~. Only two lines 
were observed which were identified with the 
two prominent lines C and D of the other spectra. 
The spacing between these two lines, as given in 
Table I, is about 4.6 times the corresponding 
spacing of the other acceptors studied. Strictly 
on the basis of effective-mass formalism, Zn’, 
having a doubly charged center, is expected to 
have spacings between the excited states which 
are 4 times larger than those for a singly charged 
center. However, the excited states concerned 
in the cases of the C and D lines of Zn” have 
binding energies already comparable to the bind- 
ing energies of the ground states of group III im- 
purities. A somewhat larger spacing between 
the C and D lines of Zn~ is not unexpected. In 
the case of Cu™, we have not yet seen the ex- 
citation lines. However, indications of the ioni- 
zation absorption have been observed. 

In the case of group III impurities, we have 
estimated the ionization energy by adding to the 
energy of the A line (ground state to 2p* transi- 
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tion) the theoretically calculated binding energy 
of the 2p* state, 0.0006 ev.° This procedure 
minimizes possible error in the theoretical esti- 
mate as the 2p* state is the highest excited state 
calculated. The A line has not yet been observed 
for Cu and Zn. The ionization energies for 
these impurities may be estimated by adding to 
the energy of the D line the binding energy of the 
2p” state which was obtained from the spectra 

of the group III impurities. Since the C-D spac- 
ing is 4.6 times larger for Zn’, we may estimate 
the ionization energy by adding to the D line a 
proportionately larger binding energy. The value 
obtained is close to the reported thermal ioniza- 
tion energy, but must be regarded as only an ap- 
proximate estimate. 

The effect of the magnetic field on the excitation 
lines has been investigated. The observations 
confirm the assignment of the lines. The dashed 
curves in Figs. 1 and 2 show the results obtained 
at a field of 16 kgauss. The D line splits into 
two components in all cases, Cu, Zn, and Zn’. 
The splitting is proportional to the magnetic field 
and is comparable to that observed for the D line 
of boron impurity. For Cu, Zn, and B, the high- 
energy component appears to be broader, and 
there is indication that this component may further 
split. Studies with higher fields and polarized 
light are in preparation. Such measurements re- 
ported recently for aluminum in silicon® show 
rather complex Zeeman effects. The situation 
in silicon is complicated by the fact that the 
separation of the third valence band by the spin- 
orbit interaction is very small. 

The C line in the Cu and Zn spectra broadens 
and shifts to higher energy. It resolves into two 
peaks at higher magnetic fields. The separation 
of each peak from the zero-field position of the 
line appears to increase quadratically with the 
field. The effect appears to be identical with that 
of the C line of boron. The Zeeman effect of the 
C and D lines reinforces the conclusion that the 
excited states of Cu and Zn are similar to those 
of the group III impurities, and indicates that the 
observed effect comes mainly from the excited 
states. In the case of Zn’, the shift of the C line 
is much smaller as can be seen by comparing 
Figs. 1 and 2. This is as expected. With a doubly 
charged center, the dimensions of the wave func- 
tions of Zn™ should be about one-half of that for 
singly charged centers. The quadratic Zeeman 
effect being a measure of the area of the wave 
functions, it might be expected to be about four 
times smaller for Zn~. This is consistent with 
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the observation. 

It is interesting to note that the absorption lines 
of impurities are very much sharper in german- 
jum than in silicon. The full width at half absorp- 
tion is about 0.19 x 107° ev for the observed lines 
of Cu and the group III impurities. The actual 
width is probably even smaller, since the ob- 


served width is comparable to the resolution used. 


This is to be compared with the width of 1x10 
ev reported for boron lines in silicon.” The 
broadening of the absorption lines has been 
treated’»® theoretically on the basis of the scat- 
tering of the bound hole by the acoustical vibra- 
tion of the lattice. According to these treatments 
which used the hydrogenic approximation for the 
impurity, the width is determined mainly by the 
broadening of the ground state, and for german- 
ium and silicon the ratio of the widths is expected 
to be of the order of the reciprocal of the Bohr 
radius squared. On this basis, we would expect 
the lines of the group III impurities to be 6 to 7 
times narrower in germanium than in silicon. 
The observation shows that the ratio of the line- 
widths in these two materials is at least as large 
as this estimate. The lines of Zn™ are consider- 
ably broader, about 0.43 x107° ev in width for 

the D line. The width is considerably larger 

than the instrumental resolution. This result 
further substantiates the expectation that the 
broadening of the ground state increases with de- 
creasing orbit dimension. In this connection, the 
failure to observe the excitation lines of Cu~ 
might be caused by large broadening. The orbit 
of the hole in the ground state may be quite small 


for Cu as indicated by the large ionization energy 
of 0.33 ev. 

One additional point of interest is the difference 
in the relative strengths of the various lines for 
different impurities. This is particularly pro- 
nounced in comparing the spectrum of Zn with 
that of the group III impurities. The C line is 
much weaker than the D line in the case of Zn’. 
Noticeable variations can be seen even among 
Cu, Zn, and the group III impurities, which 
should reflect mainly the difference in the ground 
states. Accurate data have yet to be obtained. 
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HYPERFINE STRUCTURE OF THE F CENTER IN LiF* 


W. C. Holton, f H. Blum, f and C. P. Slichter |! 
Department of Physics, University of Illinois, Urbana, Illinois 
(Received August 1, 1960) 


The F center in LiF is one of the color centers 
in alkali halides which has been most throughly 
studied experimentally by magnetic resonance’~® 
and theoretically.* To date there is serious dis- 
agreement between experiment and theory. Lord? 
originally reported resolved structure of the elec- 
tron spin resonance, which he interpeted as 
arising from the hyperfine coupling of the electron 
to the six Li’ atoms nearest to the vacancy. Al- 
though he failed to see all of the 19 lines predicted 
by this model, he later supported his contention 


with the direct measurement of the hyperfine con- 
stants by the electron-nuclear double resonance 
technique (ENDOR).*° Stimulated by this work, 
Gourary and Adrian® calculated the hyperfine 
coupling with the nearby nuclei. While the agree- 
ment with the nearest lithium atoms was satis- 
factory, the prediction for the nearest fluorine 
atoms was too large by a factor of 30. The pic- 
ture was further complicated when Kim, Kaplan, 
and Bray‘ reported that they could resolve at 
least 10 more lines than the 19 predicted. Hyde® 
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has even suggested that the resolved structure is 
not due to the F center. 

We have investigated the paramagnetic reso- 
nance of x-rayed LiF crystals containing F and 
Mcenters. This Letter reports the observation 
of 35 lines in electron spin resonance (ESR), 
together with new electron-nuclear double reso- 
nance data which enable us to predict the ob- 
served electron spin resonance in great detail 
and to resolve the discrepancy between theory and 
experiment. 

The samples were obtained from Harshaw, and 
either used without treatment, or annealed in an 
inert atmosphere, or left exposed to the atmos- 
phere for over a year. They were then exposed 
to 50-kv x rays at room temperature for varying 
lengths of time ranging from 1/2 hour to 15 hours. 
Optical density (OD) measurements show that the 
F-center concentration is varied from 10'7/cm° 
to greater than 10'*/cm! while the ratio OD;,/ODy 
is varied from 600 to ~3.5. The electron spin 
resonance is observed both at room temperature 
and at 1.3°K. The crystal is located midway on 
the narrow side wall of a TE 101 mode microwave 
cavity which is slit vertically on the center line 
of the broad face, the two halves electrically 
insulated to allow rf penetration from the rf coil 
wound on the exterior of the cavity for ENDOR 
experiments at 1.3°K. The X-band paramagnetic 
resonance spectrometer consists of a transmitter 
klystron employing a modified Pound i.f. frequency 
stabilization scheme and a high-sensitivity super- 
heterodyne receiver. 


With H, parallel to the [100] direction, the elec- 
tron spin resonances of all specimens exhibit 35 
lines with a spacing of ~14 gauss, uniformly 
decreasing in amplitude from the center of the 
line and symmetric about g=2.0006+ 0.0006. A 
typical slow-passage derivative record is shown 
in Fig. 1. To determine whether or not two reso- 
nances are present, one broad and structureless, 
the other possessing structure, we have checked 
the ratio of the resolved structure amplitude to 
the broad resonance amplitude of each sample. 
It is constant to within a factor of 2, thus essen- 
tially independent of the F/M concentration ratio 
for the range of concentrations we have used. 
The resolution is improved tenfold when H, is 
parallel to [111] and decreased twofold when H, 
is parallel to [110]. The over-all shape of the 
integrated experimental curve is nearly Gaussian 
with a width of 85+ 10 gauss from peak to center 
of the derivative. 

An attempt to explain the excessive number of 
lines by considering either forbidden transitions 
or an isotope effect is not successful when one 
uses Lord’s values of the hyperfine couplings. 

Electron-nuclear double resonance signals can 
be grouped, each group arising from a shell of 
nuclear sites which are crystallographically 
equivalent with respect to the F-center negative 
ion vacancy and whose distance from the vacancy 
is VN times the lattice parameter, where N 
is an integer. For a given shell N, there are two 
sets of ENDOR lines arising from cases where 
the field of the electron at the nucleus either aids 
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or opposes the external field. (We call the two 
cases respectively “up” and “down.”) To check 
that lines were properly identified as “up” or 
“down, ” we have studied the variation of their 
frequency with external field. We verified that 

the Shell 1 lines were correctly identified. A 
search for high-frequency ENDOR lines not pre- 
viously reported, however, reveals a complete 
set of lines having the symmetry of Shell 2, the 
nearest fluorines. The nuclei giving rise to these 
lines are identified as fluorines from the depend- 
ence of the ENDOR frequency on H, and by the 
separation of the “ups” and “downs” by twice the 
Larmor frequency in the static field (see Fig. 2). 
From the magnitude of the hyperfine splitting we 
conclude these are in fact Shell 2 and that Lord’s 
Shell 2 lines are Shell 8. The ENDOR pattern is 
usually computed using first-order perturbation 
theory. If the hyperfine coupling is treated to sec- 
ond order, subsidiary splittings arise." Such a set 
is illustrated by the 4 lines near 31 Mc/sec which 
arise from the Shell 2 fluorines. The calculation 
predicts a separation of 0.5 Mc/sec with the inten- 
sity ratios 1:4:4:1 as observed. In addition, the 
ENDOR lines were observed to obey the correct 
angular dependence of Shell 2 when the crystal is 
rotated about a [110] axis as well as a [100] axis 
perpendicular to H,. 

ENDOR lines are identified through Shell 8 
(there is no Shell 7 by this scheme of numbering 
shells); the hyperfine constants a and 6 are given 
in Table I together with the theoretically predicted 
values of Gourary and Adrian. The Shell 1 and 
Shell 2 values are in rough agreement with the 
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Table I. Hyperfine constants in megacycles/sec for 
the first seven shells of nuclei surrounding the F center 
in LiF. 








Shell Nucleus a b 

1 Li’ 39.06 + 0. 02 3.20+0.01 
(theo. = 50) (theo. =2.8) 

2 F’® 104.9440.02 14.96+0.02 
(theo. =61) 

3 Li’ 0.50 + 0.01 0.68 + 0.02 

4 F’® 0.48+0.01 1.12+0.01 

5 Li’ 0.27+0.01 0.28+0.01 

6 F’® 0.88+0.01 0.69 +0.01 

~ y* 1.34+0.01 0.56 +0.01 





predicted values. 

Using the measured Shell 1 and Shell 2 hyper- 
fine values and assuming a rectangular broaden- 
ing function to take account of the other shells, 
we have calculated the detailed ESR line shape 
for H, parallel to [111]. The result is shown in 
Fig. 3 together with an integrated experimental 
record and a Gaussian. The calculation predicts 
an average resolved structure separation 6% 
greater than the Shell 1 a value. This effect is 
clearly seen in the integrated experimental ESR 
curve. Both the experimental and calculated 
curve are nearly Gaussian. The root second 
moment, calculated from the measured hyper- 
fine constants, is 75.7 gauss; the measured root 
second moment (from derivative peaks assuming 























FIG. 2. ENDOR lines 
showing portions of Shells 
1 and 2. 


NA 
AL TT 
MN 
a 


sa HP 


nh 
vn 
wt 
wu 
mt 





en 








199 











VOLUME 5, NUMBER 5 


PHYSICAL REVIEW LETTERS 








SEPTEMBER 1, 1960 











T i . . = © 


| ie. t T | 


eee SS ee Se ee 


He along [111] af 


Q, ;7* 13.95 gauss 


= ee aD 





-—-—- Gaussian 





I clhnsilbiil 


Theoretical Calc 


Experimental 


FIG. 3. Comparison of the 
integral of an experimental ESR 
curve for Hy parallel to [111] 
with the calculated ESR line 
shape and with a Gaussian curve, 
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a Gaussian) is 85+ 10 gauss. We conclude that 
substantially the entire resonance arises from 
F centers, that the observed ENDOR lines ac- 
count completely for the experimental line shape 
and are in substantial agreement with the calcu- 
lations of Gourary and Adrian. 

Our new increased second moment for the F 
center is larger than that attributed by Lord both 
to the F center alone and to the combination of 
the F and Mcenters.” If we are correct, Lord’s 
detection of the existence of the M-center reso- 
nance by second moment methods is inconclusive. 

We would like to thank Professor Robert Maurer 
for making available his Cary spectrophotometer 
for determination of the optical spectra, and Dr. 
J. Kingsley for his helpful advice in its use. A 
substantial amount of the numerical analysis of 
the double-resonance data was performed by Mr. 
Nils Fernelius. The numerical computation of 
the theoretical line shape was performed by Mr. 
James Bushnell. We wish to thank Dr. George S. 
Newell who built the microwave spectrometer in 
conjunction with two of the authors (WCH and HB). 
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DIFFUSION OF LITHIUM THROUGH TUNGSTEN 


G. M. McCracken and H. M. Love 
Department of Physics, Queen’s University, Kingston, Ontario, Canada 
(Received August 2, 1960) 


The diffusion coefficients of Li’ and Li® in 
tungsten have been measured. The lithium is 
deposited on one surface of a tungsten ribbon and 
the atoms which diffuse through the tungsten to 
the opposite face are ionized and identified in a 
mass spectrometer. 

In a typical experiment a circular lithium de- 
posit about 100 atoms thick and 0.3-mm diameter 
was deposited from a collimated positive-ion 
beam from a thermal ion source containing the 
lithium isotopes in the natural proportions. The 
tungsten ribbon was 1.0 mm wide and 0.030 mm 
thick. After the deposition the ribbon was elec- 
trically heated to a specific temperature in the 
range 1000°K to 1400°K. The high temperature 
caused diffusion of the atoms to the other face of 
the tungsten ribbon which served as the positive- 
ion source surface for a 60°, 25-cm radius, 
direction-focussing mass spectrometer. The 
efficiency of ionization of the atoms reaching the 
second surface of the ribbon was high, as the 
work function of the surface was greater than the 
ionization potential of the diffusing atoms. When 
the ion current was observed as a function of 
time, a characteristic curve was obtained as 
shown in Fig. 1(a). 

A solution to the diffusion equation 8C /at =DV*C 
can be found for this particular arrangement if 
the diffusion through the filament of thickness L 
is characterized by a single diffusion coefficient 
D independent of concentration. It is assumed 
that the initial deposit is circular of radius b, 
and of thickness much less than L, and that every 
atom reaching the second face of the ribbon evap- 
orates from it. The current, J, integrated over 
the second face at any time ¢ after the diffusion 
process begins is given by 


I=ALD™"}*t-% exp(-L?/4D?), 


where A is a measure of the initial concentra- 
tion. From this equation it-can be seen that the 
time to reach the maximum current is L?/6D. 

D can therefore be found from this time. The 
value of D can also be found from the slope of 

the curve obtained when log(/¢*”) is plotted against 
1/t. Figure 1(b) shows such a plot for the exper- 
imental observations of Fig. 1(a). It can be seen 
that in the region around the peak current the 
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FIG. 1. (a) The current of Li’ diffusing through a 
polycrystalline tungsten sample at 1160°K as a func- 
tion of time. From the time to the peak, the diffu- 
sion coefficient, D, is obtained as 2.4x10~* cm? sec”. 
(b) The same results are now plotted in the form 
log(/t*?) vs 1/t for comparison with the theoretical 
straight line. The slope gives a value of 2.2x10~° 
cm? sec™ for the diffusion coefficient. 


experimental points lie close to a straight line. 
The values for D obtained by these two methods 
agree closely. 
Measurements of the diffusion coefficient as 
a function of temperature gives the usual tem- 
perature dependence D =D, exp(-¢€/kT), as seen 
in Fig. 2, where e€, the activation energy, is 
2.42 + 0.05 ev and D, for Li’ is 82+ 40 cm? sec™. 
The values of the diffusion coefficient and of 
the activation energy indicate that lattice diffu- 
sion through the ribbon is being measured rather 
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FIG. 2. The diffusion coefficient of Li’ in poly- 
crystalline tungsten as a function of temperature. 
The slope gives a value of 2.42 ev for the activation 


energy. 


than grain boundary or surface diffusion. This 
supposition was confirmed by first varying the 
ribbon width and then by using a single-crystal 
ribbon.’ Results for the diffusion coefficient in 
both cases were of the same order of magnitude 
as the original values and within the error in D 
due to uncertainties in the temperature and thick- 









ness of the ribbon. This indicates that lattice 
diffusion is playing the dominant role. 

The discrepancy between the observed and 
predicted currents at long times is to be expected 
as the evaporation from the first face is not taken 
into account in the analysis. The high current at 
short times, which is observed both with single 
crystals and with polycrystalline ribbon, could 
be due either to a surface diffusion contribution 
or to the diffusion coefficient being dependent on 
concentration. ' 

This method of measuring diffusion coefficients 
allows a very simple comparison of the diffusion 
rates of Li® and Li’. By switching back and forth 
between the two isotopes during a run, both curves 
can be obtained by interpolation. An average 
value of the ratio of the diffusion coefficients ob- 
tained was D,/D, =1.07+0.02 which is in good 
agreement with the ratio of the inverse square 
roots of their masses, (M,/M,)”*=1.080. Schoen’ 
has shown that the ratio of the diffusion coeffi- 
cients Dy and Dg of two isotopes of masses M, 
and Ma gives an indication of the type of lattice 
diffusion taking place. The value of Dg /Dg 
=(Mg/Mq)” is consistent with the interstitial 
diffusion which is to be expected from the sizes 
of lithium and tungsten atoms and the tungsten 
lattice spacing. 

We are indebted to Dr. P. A. Puhach of this 
Department for the theoretical analysis, and to 
the National Research Council of Canada for 
financial assistance. 
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ELECTRIC - MONOPOLE-ENHANCED 2* ° ~2* TRANSITIONS IN Th”? AND U*8*T 


F. E. Durham, ”* D. H. Rester,* and C. M. Class 
.The Rice University, Houston, Texas 





(Received June 27, 1960; revised manuscript received August 1, 1960) 


In the course of a survey for internal conver- 
sion electrons from the Coulomb-excited isotopes 
Th?*°s?52 and U*™»29*,238 a number of transitions 
were observed which are interpreted as electric 
monopole (£0) enhanced transitions between 2* 
beta-vibrational states and 2* first excited rota- 


tional states. The occurrence of E0 enhancement 
in these cases is in itself important for the iden- 
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tification of the parent Coulomb-excited states as 
beta-vibrational, since the K-selection rule for- 

bids monopole-enhanced transitions from y-vibra- | 
tional states. However, to use these resultsto _ 
test current theories of monopole enhancement** 
requires that the £0 transition probabilities and 
the ratio of the E0 to E2 transition probabilities 
for these states be determined. To do this it is 
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necessary to know not only the partial reduced as 2*’ ~2*, originating in a 775-kev §-vibrational 
transition probabilities e(e~)B(E2)4 given by the state and a 788-kev y-vibrational state, respec- 
present measurements, but also the partial re- tively.° The magnitude of the barely detectable 
duced transition probabilities €(y)B(E2)g. Since electron yield from the 788-kev transition was 
the latter results are available at present only estimated also, making it possible to determine 
for the states in Th**? and U***, the following dis- B(E2; 788 ~ 50)/B(E2; 788-0). With this ratio 
cussion is confined to these cases. and the assumption of E2 character for the 738- 


The internal conversion spectra obtained from kev transition, the value of B(E2; 788 ~0)=3.8 
Th**? and U*** bombarded with 5- Mev protens are x B(E2)sp was derived without recourse to any y- 
shown in Fig. 1, and the results derived from the ray data.® To obtain the total B(E2) value for the 
measurements are given in Table I. For U*** the 775-kev transition, however, requires use of the 
991-kev transition is assigned as the EO-enhanced y-ray data. The y rays from Th*™ were first 
g*’~2* cascade from the 1036-kev $-vibrational analyzed as if originating from a single state at 
state. The E2 crossover transition to the ground 790 kev, and B(E2)q=4.3B(E2)sp was obtained.’ 


state is not seen. The partial B(E2)g value de- Monopole conversion was not taken into account. 
rived from these measurements added to that ob- This result must now be interpreted as due to the 
tained from the y-ray measurements‘ leads to the sum of y-ray contributions from both the 775- and 
tabulated total B(E2)q=2.6B(E2)sp. Two transi- 788-kev states. To split off the part due only to 
tions are seen in the Th*** conversion spectrum, the 775-kev state, B(E2; 788-0) as determined 

a strong monopole-enhanced 725-kev transition from the conversion electron measurements was 


and a weak 738-kev transition. Both are assigned subtracted from the above value leaving a net 
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FIG. 1. The internal conversion electron spectra from Th”? and U* targets bombarded 
with 5-Mev protons are shown. Self-supporting foils of natural Th (5 mg/cm?) and natural 
U (57 mg/cm?) were used. The inset shows a rerun of the region of the K-conversion lines 
of the weak 738-kev transition in Th™?, with improved counting statistics. 
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Table I. Experimental results and derived quantities relating to E0 enhancement in Th”? and U™*. 








Isotope Th? uz 
Method A Method B 
E,*’ (kev) 775 788 775 788 1036 
e(e)BE2) , re is ‘ i 
‘contientinkele . ’ , eee 0.20 
BED 0.90 0.038 0.90 
sp 
€(y)B(E2) 
d b c c 
oe 0.5 3.8 2.2 2.1 2.4 ' 
B(E2),, 
+’ + 
BiE2;2 — 2) 10% +4 104 104 
B(E2; 2° — 0) 7 3.523 7 $.621.3 7 
BE2) ; 
BED, 1.4 3.8 2. 8(3.1) » 2.6 
p 
a 4.9 E2 1.1 E2 0.16! 
+ +! + 
u(2 ;2 ~-2) 4.9 o& 1.1 o& 0.15 
Wi(E0,e";2° ~2°) 2.010"! 0 1.0x10"! 0 5.8x10"° 
2 3.4x10!? eee 3.4x10!? see 6.1510"? 
p 0.25 0 0.17 0 0.097 





@Here €(e~) refers to conversion in the 2*’—~ 2* transition only. 

brhis value is obtained by subtracting e(y)B(E2; 788— 0) from the value, 4.3, reported by Stelson and McGowan 
for the sum of the y-ray partial B(E2) values. 

°For Th”?, ¢(y) includes all the y rays from the 2*’ state. 

AThis ratio is the theoretical value. 

©The value 2.8 was obtained by Stelson and McGowan who used the electron data of E. M. Bernstein (private 
communication); use of the present data gives 3.1, which was used to determine p(2*’; 2*’— 2+), 

fhe y rays from the 2*’— 4* transition in U™* were not included in e(y)B(E2), and hence were not included in 
the derivation of a+’ ,+ for U™*. 

8E2 character assumed. 


for €(y)B(E2; 775-0) =0. 5B(E2) gp- This value which they assume to be applicable.’ By this 
together with €(e-)B(E2)q=0.9B(E2) sp yields means they obtain ¢(y)B(E2), for the 775- and 
B(E2; 775~0)=1.4B(E2)gp. These results are 788-kev states separately, and finally B(E2; 788-0) 
summarized in Table I under the heading“Method and B(E2; 788 —50)/B(E2; 788-0). This informa- 
A,” and represent the maximum utilization of tion is summarized under the heading “Method B” 
the conversion electron data. in Table I. 

More recently Stelson and McGowan‘ also have These results have been used to evaluate the 
taken into account the presence of the 775-kev quantities of theoretical interest. The internal 
state in analyzing the y-ray spectrum of Th*™. conversion coefficient a for the E0-enhanced 


While the 725-, 738-, 775-, and 788-kev y rays 2+’ . 2* transitions is obtained by noting that 
are not resolved in their measurements, they 

have decomposed the composite peak by making 
use of the known intensity of the 612-kev transi- 








tion to the 4* state from the 775-kev state, to- picasa Oe 1. WE; 2*’ 0) 

gether with the theoretical ratios <(y)BE2) WE2, 7:2" 7m 9*) 

B(E2; 2° —0):B(E2; 2° —2°):B(E2; 2° 4") WE2,7;2" ~4°) s e) 
= 1; 40; WB, (1) WE2,y; a*’ 2°) 








8 = 0) 
ma- 
iB” 


(2) 
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The transition probabilities W(E2,y) are known 
only for the 2+’ 0 transitions; therefore, it is 
necessary to use the theoretical ratio (1) to deter- 
mine the relative transition probabilities for use 
in (2). The relation 


+’ +/ —+ W(EO e732" <3") 
y(2* 32° —2*) 2S gE 
W(E2,y;2" ~2°) 





=a - Qs, 


where a, is the theoretical E2 conversion coef- 
ficient, gives the desired ratio of transition prob- 
abilities. The monopole transition probability 
W(E0,e~; 2+’ ~2*) was determined in turn from 

(1) and the experimental B(E2; 9” « 0) values. 
The nuclear monopole strength parameter, p, is 
defined by 


W(E0, e~) =Qp?, 


where 2, the electronic factor, is given graphi- 
cally.” 

Reiner*® has treated the problem of monopole 
enhancement for deformed nuclei and has derived 
aformula for the ratio of the E0 to E2 transition 
probabilities from 8-vibrational states. For Th**? 
his formula leads to .(2*’; 2+’ +2*)=3.8, a value 
which is intermediate between the experimental 
values obtained by Methods A and B. For U*", 
however, Reiner’s formula gives .(2*’; 2+’ ~2*) 
=1.4 which is a factor of ten larger than the value 
obtained experimentally. The explanation of the 


discrepancy is not apparent in terms of the simple 
vibrational model used by Reiner, since Th**? 
and U*** have similar atomic numbers, equilibri- 
um deformations, and vibrational energies which 
are the only parameters entering into the calcu- 
lation. The values obtained for the nuclear 
strength parameter p are close to the approxi- 
mate value of 4 predicted by Church and Weneser 
for 0* ~0* transitions from vibrational states in 
spheroidal nuclei. Again, however, it is not 
clear why the value of p should differ by a factor 
of two between Th*** and U*5*. Perhaps one sees 
in these discrepancies the influence of the dif- 
ferent ground-state configurations. 





tThis work was supported in part by the U. S. Atomic 
Energy Commission. 

*Now at the Physics Department, Tulane University, 
New Orleans, Louisiana. 

‘kK, Alder, A. Bohr, T. Huus, B. Mottelson, and 
A. Winther, Revs. Modern Phys. 28, 432 (1956). 

*E. L. Church and J, Weneser, Phys. Rev. 103, 
1035 (1956). 

5A, S. Reiner, thesis, University of Amsterdam, 
1958 (unpublished). 

‘Pp, H. Stelson and F. K. McGowan (private com- 
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NONLOCAL POTENTIAL AND ALPHA-DECAY BARRIER PENETRABILITY 


Manujendra Lal Chaudhury 
Department of Physics and Mathematics, Indian School of Mines and Applied Geology, Dhanbad, India 
(Received July 12, 1960) 


Since the early works on the alpha disintegra- 
tion problem, the barrier for the alpha emission 
has traditionally been taken to be purely Coulomb- 
ian with an abrupt cutoff at the nuclear boundary 
defined by some appropriate A™’ law. Alpha 
decay was first considered’ by the present author 
by assuming that the potential barrier is not pure- 
ly Coulombian and should be taken as 2(Z -2)e?/r 
-V, where V is defined as the short-range inter - 
action between the emitted alpha particle and the 
residual nucleus. Subsequently, the form of the 
interaction, V, which is superposed on the 


Coulomb field, has been taken differently by dif- 
ferent authors? ° in their calculations of the bar- 
rier penetrability. 

It is, however, worth while to note that in the 
above calculations of the penetrability factor, 
the nuclear part of the barrier has been taken to 
be purely static, while the careful analyses of 
the recent experimental findings strongly sug- 
gest®~® the velocity dependence of both the nucleon- 
nucleus and the alpha-nucleus interactions. On 
the other hand, owing to the uncertainties in the 
knowledge of the effect of the intranuclear dynamics 
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on the half-lives of the alpha-emitting elements, it 
is important that the calculation of the penetrabil- 
ity factor be made with reasonable certainty by 
taking all such points into consideration. 

The purpose of the present work is therefore to 
study the effects of the velocity-dependent char- 
acter of the alpha-nucleus interaction on the mag- 
nitude of the alpha-decay barrier penetrability 
factor. For this purpose we have taken an inter- 
action kernel possessing a nonlocal part repre- 
sented by a 6 function approximated by a Gaus- 
sian exponential, so that 


V(t, r’)= -vos( 


_ 7792p -3 iy r-r' 2 


where V, and the form function, f (7), of the static 
part are taken from Igo’s potential, viz., V,=1100 
Mev and f (7) =exp[-(7 - 1.17A™)/0.574], r being 

in fermis; and the extent of the nonlocality, 5b, 

is given by Frahn’® for nucleons as 5 =0.902 x107** 


P+r’ 


: (1) 











)o,e -F’), 
and 


(1.1) 














cm, which is also assumed by us to be the range 
of nonlocality of nuclear potential for a@ particles, 
Now using the above potential in the integro- 
differential equation analogous to that taken by 
Frahn and Lemmer," we find the radial part of 
the Schrédinger equation for the emitted alpha 
particle to be 





Bu /dr* + (2u /h)\[e(r)E - (? /2u)U(1+1)/r7 
-2(Z - 2)e7e(r)/r+Vof (r)e(r) ey 
+ne(r){f "(r)/4 -f'/2r+f'u,'/u} =0, (2) 


where 


€(r)=[1+nf(r)]"* and n=(ub?/2h?)v,. (3) 


Now in connection with the studies of nuclear 
energies, Green’” has shown that the terms cor- 
responding to the part within second braces in (2) 
contribute negligible shift of the eigenvalues. 
Therefore, neglecting in Eq. (2) the terms in- 
volving derivatives of f(r), and using the WKB 
method of solution for simplicity, we find that the 
penetrability factor for the ground state transi- 
tions is given by 





Yr 


To _2)\p2 
exo} -2( 5) f [22 abe S07 vapirreir)-Betr] ah (4) 


v. 
1 


The integral in the exponential has been cal- 
culated numerically by using Simpson’s rule and 
taking 120 strips. The turning points 7; and 7, 
remain unaffected as for a static potential since 
the modifying factor ¢€(v) in the integral is com- 
mon to all the terms. 

It is interesting to note that the effect of the 
nonlocal interaction is quite remarkable. For 
example, in the case of ,,,Fm*™, calculation 
gives a value for the penetrability factor of 
5.848 x10~%*, which is about double the value ob- 
tained with the corresponding static potential, 
viz., 3.524x1074. 

Numerical calculations for the whole series of 
the alpha-active elements and for their excited- 
state transitions are in progress and a detailed 
report will be published in due course. 
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CHARGED PARTICLE AND TOTAL REACTION CROSS SECTIONS FOR PROTONS AT 9.85 Mev* 


Verena Meyer? and Norton M. Hintzt 
School of Physics, University of Minnesota, Minneapolis, Minnesota 
(Received June 27, 1960) 


The optical model can predict elastic scatter- 
ing, polarization, and total reaction cross sec- 
tions. Recent parameter studies’ indicate that at 
least a necessary condition for obtaining a unique 
set of parameters is the existence of data cover- 
ing all three quantities at the same energy. Un- 
til now there has been no energy (or bombarding 
particle) for which such complete data are avail 
able. 

Near 10 Mev, both differential elastic?» * and 
polarization measurements* have been made for 
many nuclei with protons, and thus reaction 
cross sections are particularly desirable. In 
addition, there are reasons to believe? that it is 
just at charged particle energies in the vicinity 
of the Coulomb barrier that reaction cross sec- 
tions will be most useful in resolving the ambi- 
guities in the optical model parameters. This 
note reports measurements of the sum of the 
cross sections for all reactions in which a 
charged particle is emitted. This quantity, 
which we designate o(,q), when added to the 
total (p,m) cross section gives essentially the 
total reaction cross section, Op. The values of 
0p which we obtain are shown to disagree with 
current optical model predictions. The signifi- 
cance of this disagreement is discussed. 

Direct attenuation measurements have been 
made of op with protons of 34 and 61 Mev.*»° 
These methods become progressively more dif- 
ficult as the proton energy is reduced. At 10 
Mev, however, it is feasible to measure all of 
the significant partial reaction cross sections, 
ie., o(p,n), o(p,p’), o(p,a), etc., and sum to 
obtain op. Reactions in which two or more par- 
ticles are emitted are either energetically for- 
bidden or greatly suppressed because of the 
Coulomb barrier. The experiment reported here 
was designed to measure the energy and angular 
distribution of all charged particles emitted from 
the nucleus down to the lowest energy allowed by 
the Coulomb barrier. Separate measurements 
are then needed of the (p,m) cross section to 
obtain Op- The experiment was performed using 
9.89-Mev protons from the Minnesota linac. A 
CsI detector and pulse-height analysis were used 
to obtain complete charged particle spectra at 
10° intervals from 20° to 150° for a number of 


elements. Figure 1 shows a spectrum for Zn® 

at 150°. The shaded region shows the background, 
taken with an absorber over the detector to stop 
all charged particles. 

The spectra were integrated over energy (leav- 
ing out the elastic peak) and angle to obtain o( p,q) 
for each element. The results are given in 
Table I. The angular distributions of the non- 
elastic charged particles were found to be iso- 
tropic except for a slight forward peaking due to 
a small direct component in the inelastic protons 
to the first excited state of the target nucleus. A 
complete presentation of the spectra together 
with an analysis of the energy and angular distri- 
butions will appear in a later paper. 

It can be seen from Table I that there is a con- 
siderable variation in the values of o(p,q) from 
element to element. This is to be expected from 
a statistical theory of the decay in which the ra- 
tio of inelastic scattering to neutron emission is 
strongly dependent on the (p,m) Q value and on 
the Coulomb barrier height. 

To obtain op, the values of o(p,n), o(p,y), and 
Oce (compound elastic) must be added to o(p,q). 
Measurements of o(p,y) near 10 Mev give values 
of 1-2 mb,’ so this quantity can be neglected. 
Compound elastic cross sections for protons 
have not yet been measured. An upper limit to 
See would be approximately the isotropic part 

















121 
> 

= 

se, 10h 

Q 

3 

= BL 

| 

8 6| ] 
3 

2 

a 4+ 

3 

So a } 





; + ¢ &4eett es Se 
Energy — Mev 


FIG. 1. Charged particle spectrum for Zn® at 150°. 
Energies are in the laboratory system. 
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Table I. Charged particle and total reaction cross sections at 9.85 Mev. 



























o( p,q) o(p,n) 0,, experimental ¢,, optical model 
a R R 
Target mb mb mb mb 
Cc >228 + 15? 
+ 60 c + 60 
Al 725 * $0 >16 >741 * 35 606 
Ti 362 + 25 >371 + 304 >733 + 40 
+170 +170 
Vv >165 + 10° 552 ~ 190 >717 _ 190 
Fe 498 + 33 >367 + 40! >865 + 52 
Co 366 + 27 
Ni 590 + 30 >1238 >713 + 40 576 
Cu’ 365 + 18 4807 30 eas * 32 
Cu 155 + 10 819 + 75 974 + 76 670 to 800" 
Zn 306 + 21 640 + 80° 946 + 85 
Ag 30.7 *3 
Sn 28+ 6 425 



























© Value of o( p,m) at 6.7 Mev, reference 10. 


SDoes not include Ni*®(p, 2)Cu®. 


and L. Grodzins, Phys. Rev. 119, 276 (1960). 


of the cross section for exciting the first state 
of the nucleus. For the medium weight even-even 
nuclei this was found to be ~50 mb. For the odd- 
A nuclei Ogg is expected to be considerably less.*® 
Unfortunately only a few measurements of 
o(p,n) exist near 10 Mev,°® ** and these measure- 
ments are not in good agreement with each other. 
The values of o(p,) shown in Table I are more 
or less arbitrarily weighted averages of the var- 
ious values. For the natural isotopic targets 
the o(p,m) value quoted is the weighted sum of 
the cross sections for those isotopes for which 
measurements exist. The errors indicated cover 
the disagreements between the several experi- 
ments. The column marked “op, experimental” 
is the sum of our o(p,q) value and the value 
chosen for o(p,n). Only lower limits to op can 
be given when (p,m) measurements are not avail- 
able for all of the isotopes of the target. No at- 
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hRange of values for Cu for %)= 1.20 to 7)= 1.33, references 1 and 16. 
1 Based mostly on (P,”) value from reference 11 with correction for p°/K ratio determined by A. Schwarzschild 


2 The targets were of natural isotopic composition except for Cu and Zn. 
D some of the low-energy alphas from the (f, 3a) reaction were below the detection threshold. 


4 Does not include o(p,”) for Ti» 4% 4 (20.5% relative abundance). 
©Does not include o(,p') to states in V*! at 0.32 and 0.65 Mev which could not be resolved. 
foes not include o(p,”) for Fe™ (5.8% relative abundance). 


tempt has been made to estimate o,,, and thus 
the op values must be considered lower limits 
to the total reaction cross sections as computed 
in the optical model. It is felt that the large dif- 
ference in Op» XP, for the two copper isotopes 
may not be real but may be due to the large er- 
rors in the (p,m) measurements. For compari- 
son with the optical model, the average for the 
two copper isotopes, op = 910+ 60 mb, should be 
used, along with the value of op = 946+ 85 mb for 
Zn®*. In the optical model, the value of op for 
Zn® is expected to be within a few percent of 
those for the copper isotopes. In the schematic 
theory of compound nucleus formation, in which 
the nucleus is depicted as a uniform, completely 
absorbing sphere, a radius constant of 7,=1.75 
is needed to give op = 946 mb for Zn®. This is 
a considerably larger value than has been used 
in fitting (p,”), (a,m), etc., excitation func- 
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tions,** but in these cases a significant part of 
the reaction cross section has been omitted, 
namely, the inelastic scattering.** 

The relevant optical model values of op are 
those of Glassgold et al.** using a central poten- 
tial and those of Nodvik and Saxon’ calculated 
with the inclusion of a spin-orbit term. Nodvik 
and Saxon obtain for copper, after fitting the 
elastic scattering and polarization data, a value 
of Op = 670 mb for 7,= 1.20 and op = 800 mb for 
7,=1.33. The best fit to the elastic and polariza- 
tion data is obtained for 7,=1.25, for which op 
=720 mb. The values of op obtained by Glass- 
gold, without the spin-orbit term, are very close 
to the above. In both cases calculations were 
made with the same Saxon form factor for the 
real and imaginary parts of the potential, V(r) 
and W(r). 

It is seen that the optical model values are 
from 100 mb to 200 mb lower than the experi- 
mental values for iron, copper, and zinc and in 
addition the o(p,q) values alone for Ni and Al 
are comparable to or greater than the optical 
model values for op. The inclusion of o,, in the 
experimental op would further increase these 
discrepancies. 

Preliminary calculations by Easlea’’ show that 
by putting all of the absorption at the nuclear 
surface the reaction cross section can be raised 
to 844 mb for copper with 7,=1.30, still fitting 
the elastic and polarization data. 

We feel that this value of op is still too small 
and that our data constitute strong evidence for 
the extension of W(r) beyond V(r) as has been 
predicted theoretically.** More calculations as 
well as more accurate (p,m) cross-section meas- 
urements are needed to settle this question. 

Finally it should be noted that either a surface 
peaking or an extension of W(r) to a larger ra- 
dius will lower the central value of W to one 
more in line with predictions of the Brueckner 


theory for infinite nuclear matter.” 





*Supported in part by the U. S. Atomic Energy Com- 
mission. 

tNow at Physikalisches Institut, Universitat, Ziirich, 
Switzerland. . 

tNational Science Foundation Senior Postdoctoral 
Fellow on leave at the Institute for Theoretical Physics, 
University of Copenhagen, Copenhagen, Denmark. 

‘See, for example, J. S. Nodvik and D. S. Saxon, 
Phys. Rev. 117, 1539 (1960). 

Norton M, Hintz, Phys. Rev. 106, 1201 (1957). 

3G. W. Greenlees, L. G. Kuo, and M. Petravié, 
Proc. Roy. Soc. (London) A243, 237 (1958); see also 
W. M. Gibson, D. J. Prowse, and J. Rotblat, Proc. 
Roy. Soc. (London) A243, 237 (1958). 

‘L. Rosen, Proceedings of the International Con- 
ference on the Nuclear Optical Model, Florida State 
University Studies, No. 32 (The Florida State University, 
Tallahassee, Florida, 1959). 

5T, J. Gooding, Nuclear Phys. 12, 241 (1959). 

Sv. Meyer, R. M. Eisberg, and'R. A. Carlson, Phys. 
Rev. 117, 1334 (1960). : 

™B. L. Cohen, Phys. Rev. 100, 206 (1955). 

8See, for example, G. W. Greenlees and P. M. 
Rolph, Proc. Phys. Soc. (London) 75, 201 (1960). 

°S. N. Ghoshal, Phys. Rev. 80, 939 (1950). 

103, P. Blaser, F. Boehm, P. Marmier, and D. C. 
Peaslee, Helv. Phys. Acta 24, 441, 465 (1951). 

‘tH, A, Howe, Phys. Rev. 109, 2083 (1958). 

12H, G. Blosser and T. H. Handley, Phys. Rev. 100, 
1340 (1955). 

18§. Tanaka and M. Furukawa, J. Phys. Soc. (Japan) 
14, 1269 (1959). 

‘“M. M. Shapiro, Phys. Rev. 90, 171 (1953). 

8In cases where the @ of the (p,m) reaction is very 
low so that this reaction will dominate, values of 7) as 
large as 1.7 are needed; see, for example, reference 
12. 

16a. E. Glassgold, W. B. Cheston, M. L. Stein, 
S. B. Schuldt, and G, W. Erikson, Phys. Rev. 106, 
1207 (1957). 

1'B, Easlea (private communication). 

'8See, for example, Gordon L. Shaw, Ann. Phys. 8, 
509 (1959) and L. C. Gomes, Phys. Rev. 116, 1226 
(1959). 











209 











VoLuME 5, NUMBER 5 


PHYSICAL REVIEW LETTERS 


1960 


SEPTEMBER 1, 





INVESTIGATION OF MEDIUM AND HEAVY NUCLEI BY 
QUASI-FREE PROTON-PROTON SCATTERING* 


‘ Gerhard Jacob and Th. A. J. Maris 
Instituto de Fisica and Faculdade de Filosofia, Universidade do Rio Grande do Sul, Pérto Alegre, Brasil 
(Received August 16, 1960) 


Recently, experimental groups in Uppsala,’ 
Harwell, ? and Chicago* have directly demonstrated 
the shell structure of several light nuclei by 
measuring the energy and angular correlations 
of proton pairs resulting from high-energy (150- 
440 Mev) quasi-free knockout processes. The 
binding energies and interesting features of the 
momentum distributions of the protons in the 1s 
and 1p shells and the spin-orbit splitting of the 
1p shell have been determined in this way. 

Up to now these measurements have been lim- 
ited to light nuclei, mainly because the cross 
section for knockout processes in the inner shells 
and the level distances in the upper shell decrease 
rapidly with increasing nuclear radius.‘ It is the 
purpose of this note to point out that high-energy 
quasi-free p -p scattering is not only able to 
supply new information on the light nuclei, but 
is also a promising tool for the investigation of 
medium and heavy nuclei. It is expected that the 
limitation in size of the effective interaction zone, 
resulting from the absorption of the protons by 
the nucleus, should give rise to a typical angular 
correlation pattern from which several impor- 
tant quantities may be read off directly. 

Throughout this paper the indices 0, 3, 1, 2, 

i, and/f refer, respectively, to the incoming and 
the nuclear proton before the collision, to the 
two outgoing protons, and to the initial and final 
nucleus. The symbol k refers to momentum; 

a; is the angle which k; makes | with the projec- 
tion of k, on the plane through k, and k, (see 

Fig. 1); m denotes a magnetic queiun number 
and P the parity operator. 


SPIN QUANTIZATION AXIS ~| 








MIRROR 
SYMMETRY PLANE 





FIG. 1. 
symmetry. 
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Consider as an example the coplanar (a, = 0) 
quasi-free scattering of 150-Mev protons by 
medium and heavy nuclei and choose a, and a, 
not too different from 45°. Comparing the mean 
free paths in nuclear matter of the involved pro- 
tons with the total length of their trajectories in 
the nucleus, one sees that the bulk of the protons 
which suffer only one collision originates from 
the two pole caps at the ends of the nuclear sym- 
metry axis orthogonal to the scattering plane.*»® 
The contribution of the multiply scattered proton 
pairs to the correlation cross section is smeared 
out over the whole available energy and angular 
range. It results in a smooth background in the 
measured energy spectrum on which the peaks 
caused by the relevant clean knockout collisions 
are superimposed. It is to be expected that the 
height of the aforementioned pole caps is small 
compared with their mutual distance D. The con- 
tributions of the two pole caps will, because of 
symmetry reasons, be equal in magnitude, though 
they have in general a phase difference 7. If one 
changes a, to small nonzero values (keeping a,, 
a@,, k,, k, constant), the magnitudes of the pole 
cap contributions remain approximately equal, 
but their phase difference becomes n+a@,k,D. For 
the transition to a definite final state of the nu- 
cleus, the contributions of the two pole caps will 
add coherently and the correlation cross section 
will show an interference pattern of the form 


| Gla, @,, 2, k,,k,)|? cos*[3(a,k,D+n)]. (1) 


As a function of a,, |G|? varies slowly as com- 
pared to the cos? term (which in our example has 
a period of about 10°). 

The value of 7 can be found by a symmetry 
argument. Consider the operator Qz which mir- 
rors the wave function through the plane a, =0. 
We quantize the spins in the direction orthogonal 
to this plane and note that Q, equals the parity 
operator times a rotation by an angle 7 around 
the quantization axis. The initial state is an eigen- 
state of Q, with eigenvalue P; exp[i(m;+m,)7]. 
The mirror parity operator of the final state is 
Py exp[i(my +m, +m,)7]Qy , where Q, is the 
mirror parity of the interference pattern. Parity 
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conservation results in the condition: 


Q,. = P,P exp[iAmr], (2) 


with 


Am =(m, +m») “(mn +m, +my). 
Because of fermion conservation, Am is an inte- 
ger. Thus, only the cases Qg _*2 1 occur. 

Using the single-particle model and the im- 
pulse approximation (for which, by the same 
symmetry arguments, m,+ m,-m,-m,=0), one 
obtains Qq =(-1)!+™, where / and m refer to 
the Y "ae component of the state from which the 
proton has been knocked out. 

From (1) and (2) it follows, independent of the 
experimental energy resolution, that the corre- 
lation cross section will have the form 


a(Qy, @,, @,k,,k,) cos?($ a,k,D) 


+b(a,, @,, @, ky, k,) sin*"(sa,k,D), (3) 


where a and 6 are slowly varying functions of a. 

If a measurement is made in such a way that 
many final states of different character contribute 
statistically to the cross section, a and b will 
have similar magnitudes and the a,-dependent 
interference pattern will fade out. On the basis 
of the shell model and of the nature of the reac- 
tion, one expects to be able to resolve groups of 
states which correspond to the shell-model state 
from which the proton has been ejected. But, 
because the m degeneracy cannot be resolved, 
one in general will still have a mixture of cos? 
and sin? contributions. 

For this last case, detailed calculations® have 
been made with a single-particle model, the 
impulse approximation, and a simplified distor - 
tion for the incoming and outgoing waves. After 


summing over the magnetic quantum numbers, 
the coefficients a and b are shown to be strongly 
different for small momentum transfers (nonrela- 
tivistically: a,=a,=45°; k,=k,=4Vv2k,). The 
expected a,-dependent interference pattern is 
therefore pronounced and gives by its period the 
distance between the pole caps, by its modulation 
the height of the pole caps, and by its phase the 
parity of the shell-model state from which the 
proton has been knocked out. Furthermore, it 
turns out that the orbital angular momentum of 
this state can be found from the (a,, a,) depen- 
dence of the correlation cross section. 

Finally, we want to point out two necessary 
conditions for the applicability of our general 
treatment. At a given bombarding energy the 
nucleus to be investigated must be at least so 
heavy that the heights of the contributing pole 
caps are rather small compared to the nuclear 
radius. On the other hand, these heights should 
not be too small, because otherwise the absolute 
magnitudes of the pole cap contributions become 
already very different for small values of a, and 
the correlation pattern fades out. A more quan- 
titative prediction requires distorted wave calcu- 
lation for each case. 





*Work partially supported by Comissfo Supervisio- 
nadora dos Institutos and Conselho Nacional de Pes- 
quisas, Brasil and by the U. S. National Science Foun- 
dation. 

'H, Tyrén, Peter Hillman, and Th. A, J. Maris, 
Nuclear Phys. 7, 10 (1958). Peter Hillman, H. Tyrén, 
and Th. A, J. Maris (to be published). 

?T, J. Gooding and H. J. Pugh (to be published). 

H. G. Pugh and K, F. Riley (private communication). 
3H, Tyrén and P. Isacsson (private communication). 
‘Th. A, J. Maris, Peter Hillman, and H. Tyrén, 

Nuclear Phys. 7, 1 (1958). 
5a, J. Kromminga and I. E. McCarthy, Phys. Rev. 
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PION CLOUD EFFECTS AND THE TWO-CENTER MODEL OF COSMIC-RAY JETS 


E. M. Friedlander 
Institute of Atomic Physics, Bucharest, Romania 
(Received July 12, 1960) 


Recent investigations of high-energy nuclear 
interactions'~* have shown that at energies above 
1 Gev an important part is played by peripheral 
collisions involving only one or two pions in the 
meson cloud of the colliding nucleons. In the 
present note, an attempt is made to interpret the 
two-cone structure of cosmic-ray jets along these 
lines. 

Indeed, it has been shown’:® that most jets can 
be interpreted as two bunches of particles emerg- 
ing from two centers moving in opposite directions 
each with a Lorentz factor y in the center-of-mass 
system of the colliding nucleons. Until now, y 
has been treated as a free parameter, supposed 
to be related in some way or other to the impact 
parameter and providing essentially an estimate 
of the inelasticity of the collision. Perhaps the 
most remarkable fact of the results obtained until 
now by means of the two-center model is the clus- 
ter of values of y about a mean value ~1.5 - 1.6. 

On the other hand, the pion-cloud effects de- 
tected at lower energy ought to appear even more 
pregnantly at the high energies of the cosmic-ray 
jets, owing to the very short deBroglie wave- 
lengths of the colliding particles. This means 
that at reasonable impact parameters, each of 
the two nucleons has a large chance to collide 
with one single pion in the meson cloud of its 
collision partner. The process is then essentially 
a pair of head-on pion-nucleon collisions.*»?° 
Owing to the fact that the incoming nucleon and 


its pion cloud move with the same velocity in the 
lab system (Lorentz factor y,), the kinematics of 
the two quasi-independent pion-nucleon collisions 
automatically lead to a fixed value of y at high y, 
values (say, above 10" ev). 

Let y, be the Lorentz factor for the collision of 
the incoming nucleon with the pion at rest in the 
lab system, and y, the corresponding quantity for 
the collision of the incoming pion with the target 
nucleon. We have then (see reference 10) 


¥1= (yo+ w)(1+ uw? + 2yuy,)™, (1) 


Y2= (Ly—+ 1)(1+ p?+ 2uy) , (2) 


where py is the pion rest mass (M=c=1). For 
Yo> 1 these expressions reduce to 


¥1* (yo/2u)™, (3) 
¥2% (uy,/2). (4) 


Therefore, obviously, 
11%2*¥0/2* 7," (5) 


where y, is the Lorentz factor of the c.m. system 
of the colliding nucleons with respect to the labo- 
ratory frame, 


¥:/¥2= (ot H)/(Uyo+1)=u-*, (for large y,) (6) 


Table I. Comparison of computed and observed quantities for the two-center model. All experimental y values 
are Castagnoli estimates®; for Yc2 7, they have to be corrected by a factor?» ¢ f~0.7. At lower Yc values this 
correction factor approaches unity. > Since y varies as ¥c"", the corresponding correction factor has been taken 
as ~f v2. In the range 5<y,.< 10 the interpolation formula deduced in reference 1 has been used. The expected 


values are weighted averages. 








Energy 7 7 ¥1/Y2 vi/"2 V:/"2 

range Y°) Computed Observed Corrected Computed Observed Corrected 
1<yce=5 3.7 1.37 1.35 + 0.05 5.20 5.15 +0.58 ae 
5<y-<10 7.0 1.43 1.56 + 0. 06 1.39 + 0.05 6.45 8.00 +1.18 6.40 + 0.95 
Ye? 10 44.3 1.49 1.89 +0.13 1.59 +0.11 6.67 13.70 + 2.42 9.56 + 1.69 
All jets 19.3 1.49 1.63 + 0.06 1.47 20.05 6.67 9.45 + 0.06 6.61+0.75 





8C. Castagnoli, G. Cortini, C. Franzinetti, A. Manfredini, and D. Moreno, Nuovo cimento 10, 1539 (1953). 


be. M. Friedliinder, Nuovo cimento 12, 483 (1959). 


°p. Jain, E. Lohrmann, and M. Teucher, Phys. Rev. 115, 643 (1959). 
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and 
y=(u+ I)y (1+ y2+2yy) = H(u?+u-), = (7) 


in the same approximation as above. 
Hence, at large energies we would expect 


¥1/¥2= 6.67, © (8) 


and 
y= 1.49. (9) 


At low energies (say, y.<7 or y,<100), both 
y,/Y2 and y rapidly decrease and approach unity. 

According to the above ideas, an analysis has 
been made of 78 jets (p, n, and a events)*! meas- 
ured in the Prague Laboratory. The results are 
shown in Table I. The agreement between the 
expected and the observed values and their trends 
is obvious. 

Further applications of this method may be 
thought of for the case of nucleon-nucleus and 
nucleus-nucleus collisions, where peripheral 
interactions are expected to be sometimes super- 
imposed on head-on nucleon-nucleon collisions. 


The author is very indebted to Dr. J. Pernegr, 
who kindly put at his disposal detailed data of the 
jets recorded and measured in the Prague part of 
Po-Stack No.1. 
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sions. 





PROOF OF THE MANDELSTAM REPRESENTATION IN PERTURBATION THEORY* 


Richard J. Eden 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received August 9, 1960) 


I wish to report a proof that the Mandelstam 
representation’ for a scattering amplitude is 
valid for every term in the perturbation series 
for the amplitude. The proof applies to any sys- 
tem of interacting particles that does not have 
anomalous thresholds. The latter requirement 
can be expressed by the mass conditions obtained 
in fourth order. In this Letter I will outline the 
main steps of the proof, using the equal-mass 
case as an illustration. The details are con- 
tained in two papers”»® in the course of publica- 
tion and a third* shortly to be submitted. 

The invariant energies squared will be denoted 
s, tj andu; a term in the expansion of the ampli- 
tude A(s, ¢) will be written F(s,¢). The main 
steps include the following: 

(1) The physical branch of F in a physical scat- 
tering region has the representation 





F(s, t) 
: (1 -Dao(a)y’ 77? 
smite { da,- . ‘da, 3 
e-0 4, [D (a,s, ty ~ 


where we have (I, Sec. 6),? 


D (a, s, t)=sf (a) +tg(a) - m?K(a) +ieC(a) a, 
=D(a,s, t)+ieC(a)Da.. 


Singularities (i.e., branch points) of F will occur 
when D¢ has end-point zeros or coincident zeros 
(pinching the integration contour) in each a vari- 
able as € approaches 0. End-point zeros can be 
reinterpreted by using reduced diagrams. 

(2) D(a, s, t) is negative for real positive a in 
the region s>0, ¢>0, u>0. Hence F is real in 
this region [I, Secs. 4 and 8 (F)]. 

(3) The only straight lines of singularities of 
the physical branch of F are normal thresholds, 
since they must intersect a physical scattering 
region (I, Sec. 5). 

(4) Curves of singularities I'(s,¢) of F in the 
real s,¢ plane have slope (I, Sec. 7) 


dt /ds =- f (a)/g(a), 


and they have normal thresholds as asymptotes. 
(5) From (2), (3), and (4) it can be shown (III)* 
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that there are no singularities in the region 
s<4m?, t<4m?, u<4m?, and that D is negative 
and F real in this region. It also follows that no 
curve of singularities enters the real s, ¢ plane 
in the region -4m?<t<4m?. 

(6) A single-variable dispersion relation in s 


can be proved (III) for ¢ real and with (-4m?<t<4m?’), 


by noting that for s>4m? and real a we have 
D(a, s+ie, t)#0. 


We obtain the correct physical branch of F(s, ¢) 
in the physical region s >4m? by defining it from 


F(s, t) 
n-21-1 


= limit L. =a 
e-0 “, [D(a, s +ie, t)] 





: (1 - Da [c(a)] 
ef da.-++-da 


No distortion of the real a contours is necessary 
except in the limit as € approaches 0. This proves 
that F(s,¢) has no singularities in the upper half 
s plane when ¢ is real and in (-4m?<t<4m?). It 
also defines the physical sheet of the complex 
variable s. 

(7) If a diagram, or reduced diagram, has its 
first thresholds at s=a, t=b, u=c, the limits 
+ 4m? in (5) and (6) can be extended to suitable 
combinations of a, b, andc. The discussion be- 
low can then be carried through by using these 
limits. 

(8) The curves of singularities I'(s, ¢) defined 
by ¢=¢(s), with s and ¢ real, will lie on surfaces 
of singularities (s,¢). These will be singulari- 
ties of some analytic continuation of F. Our 
problem is to prove that none of these surfaces 
correspond to singularities on the physical sheets 
of s, t, andu. This can be done by proving that 
(a) in the region s>4m?, t>4m?, all curves of 
singularities have negative slope, and (b) there 
are no disconnected complex singularities in the 
physical sheet (II).° 

(9) There are two alternative methods which 
give these results. I will illustrate the first (III) 
by examining the continuation of F(s, ¢) past t=4m? 
along any path in the upper half of the s plane and 
with ¢=t,+i¢e. A path of this type leads to the 
physical branch of F in a physical region, and it 
lies on the physical sheet. If no such path meets 
a singularity of F, our required result is true. 

If the path does meet a singularity, it may be 
either a disconnected complex singularity or a 
complex singularity coming from a spurious turn- 
ing point. By continuity, the latter must occur 





214 








at a minimum of ¢ along a curve of singularities, 
since these curves cannot disappear from the 
physical sheet except asymptotically through 
normal thresholds. If any path meets either type 
of singularity, choose the singularity on the phys- 
ical sheet for which ¢, is nearest to 4m”, if pos- 
sible. But we can then trace the singularity along 
its related singular surface © on the path s, +is, 
= s(t, +i€) to smaller values of ¢, until ¢, is less 
than 4m”. This path does not meet either t =4m’? 
or s=4m?’, since for finite ¢, the function F(s, +is,, 
t, +i€) is analytic near these thresholds. Letting 
€ tend to zero, we obtain a contradiction of our 
single-variable dispersion relation below ¢ = 4m’. 
This proves our required result that there are no 
spurious turning points and no disconnected com- 
plex singularities. 

(10) The second method (III) uses the elegant 
device of analytic completion.* The result of 
Sec. (6) is first extended to show that F(s, ¢, + ie) 
is analytic for s in the upper half plane. Then 
we define, for t, <4’, 





: 1 F(z,t, +ie)dz 
Fls,t,+i0=35 | he ee — 
Cc 


where C is a semicircle in the upper half plane. 
The contour C can now be displaced in the com- 
plex s,¢ space by continuously increasing f,. 
Provided the contour does not meet any singular- 
ities, this procedure defines an analytic continu- 
ation in the upper half s plane for all real posi- 
tive values of t,. There are no anomalous thres- 
holds, and the contour is not distorted near the 
normal thresholds since near these F(s, +is,, 

t, +i€) is analytic. Spurious turning points (min- 
ima in ¢) and disconnected complex singularities 
cannot exist, since they would lead to horns ex- 
tending into the volume swept out by C. This 
proves that the analytic completion we require is 
possible. The curves on which F(s, +ie’, ¢, +i€) 
is singular, as € and e€’ tend to 0, will have nor- 
mal thresholds as asymptotes and will have nega- 
tive slope. 

(11) A similar discussion (III) shows that 
F(s+ie,t-ie’) analytically continued from the 
real region is not singular for s>4m?, t>4m’ 
except in the limit €, «€’—0. In this limit it is 
singular only at the normal thresholds and not 
on the curves of singularities. Thus where the 
slope of I is negative, it extends to = only at 
points s+ie, ¢-ie’ that lie in nonphysical sheets. 

(12) From (10) and (11) a double application of 
Cauchy’s theorem is possible (II) and establishes 
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the validity of the Mandelstam representation for 
every term F in the perturbation series for the 
amplitude A. 

I am indebted to Mr. J. Nuttall and Dr. K. Sy- 
manzik for valuable discussions. I wish also to 
thank Dr. J. C. Polkinghorne for an outline of his 
work (in collaboration with Mr.-P. Landshoff and 
Dr. J. C. Taylor) on a different approach to a 
proof of the Mandelstam representation. 





*This work was done under the auspices of the U. S. 


Atomic Energy Commission. 

'S, Mandelstam, Phys. Rev. 112, 1344 (1958). 

*R. J. Eden, Phys. Rev. 119, 1763 (1960); denoted 
I above. 

5R, J. Eden, Phys. Rev. (to be published); denoted 
II above. 

4R. J. Eden (to be published); denoted III above. 

5G. Kallén and A. Wightman, Kgl. Danske Videnskab. 
Selskab, Mat. -fys. Skrifter 1, No. 6 (1958), and 
S. Bochner and W. T. Martin, Several Complex Var- 
iables (Princeton University Press, Princeton, New 
Jersey, 1948). 
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ABSTRACTS 


In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 








TEMPERATURE DEPENDENCE OF THE LOW- 
MOMENTUM EXCITATIONS IN A BOSE GAS OF 
HARD SPHERES. F. Mohling, Department of 
Physics, Cornell University, Ithaca, New York, 
and M. Morita, Department of Physics, Columbia 
University, New York, New York (Received 

June 24, 1960). 


The pseudopotential method is used to calculate 
the temperature dependence of the low-momentum 
excitations in a Bose gas of hard spheres to an 
order beyond that previously known. The cor- 
rection term is shown to be small compared to 
the leading term for all temperatures less than 
the critical temperature T,. The excitations 
become poorly defined for temperatures slightly 
below T,, however, since it is shown that in 
this region the excitation lifetimes are quite 
small. 


STATISTICAL MECHANICS FOR THE NON- 
IDEAL BOSE GAS. A. E. Glassgold, A. N. Kauf- 
man, and K. M. Watson, Lawrence Radiation 
Laboratory and Department of Physics, Univer- 
sity of California, Berkeley, California (Re- 
ceived May 16, 1960). 


The equilibrium and quasi-equilibrium pro- 
perties of a system of interacting bosons are 
studied from a microscopic point of view. For 
equilibrium, the model of Bogolyubov is gener- 
alized to finite temperature by using the grand 
partition function. The thermodynamic proper- 
ties and the pair-correlation function are cal- 
culated. The statistical mechanics for moving 
systems is then developed and applied to the 
problem of a rotating fluid. For quasi-equilib- 
rium, general transport equations are derived 
from first principles, independent of statistics 
and model. For the Bogolyubov model, the 
familiar two-fluid hydrodynamics is then de- 
rived, leading to the phenomena of first and 
second sound. 
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DIELECTRIC FORMULATION OF QUANTUM 
STATISTICS OF INTERACTING PARTICLES. 

F. Englert* and R. Brout, Laboratory of Atomic 
and Solid State Physics and Department of Physics, 
Cornell University, Ithaca, New York (Received 
February 25, 1960). 





Formal relations between the free energy, the 
“dielectric constant” which expresses the re- 
ponse of the system to an external perturbation, 
and the two-particle Green’s function in tempera- 
ture space are derived. Connection with pertur- 
bation expansion for free energy and for general 
admittance tensors are established. These re- 
sults are applied to a general discussion of the 
random phase approximation at finite temperature 
and it is shown that the sum on ring diagrams | 
corresponds to the calculation of the dielectric 
constant in the approximation of neglect of local 
field corrections. 

*On leave of absence from Université Libre de Brux- | 
elles, Bruxelles, Belgium. 


STREAMER MECHANISM IN FILAMENTARY 
SPARK BREAKDOWN IN AIR AND ARGON BY 
FAST PHOTOMULTIPLIER TECHNIQUES. L. B. 
Loeb, G. G. Hudson,” R. G. Westberg,! and H. C. 
Huang, Department of Physics, University of 
California, Berkeley, California (Received 
June 1, 1960). 


This study was undertaken to demonstrate by 
direct photomultiplier observation that the trans- 
ition from a positive point corona or a Townsend 
predischarge in more uniform field geometry to 
a transient filamentary arc in what is commonly 
called the electric spark proceeds as a sequence 
of luminous streamer mechanisms starting at the 
anode and crossing to the cathode at high speeds. 
Two photomultipliers with viewing slits trans- 
verse to the direction of propagation were used; 
one focused just below the anode served to trigger 
the. sweep of a 517 Tektronix Oscilloscope which 
registered the passage of luminous pulses viewed 
by the second photomultiplier at various points 
across the gap. The observations of relative 
intensity plotted against time permitted the cross- 
plotting of distribution of luminosity across the 
gap at various time intervals after the breakdown 
started at the anode. In room air at atmospheric 
pressure for a wide range of gap forms the trans- 
ition consisted in a fast primary streamer followed 
by a more intense secondary or even a tertiary 
streamer crossing from anode to cathode. On 
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arrival of the secondary or the tertiary streamer 
at the cathode, the luminous arc channel was 
established and increased in luminosity until the 
capacity was discharged. For highly divergent 
anode fields much branching of primary streamers 
occurs. The secondaries follow only a few of 
these and the crooked spark channel follows the 
path of the fastest and most effective secondary 
streamer on arrival at the cathode. These in- 
vestigations were carried over to Linde’s spec- 
troscopically pure A gas in a system with Alpert 
vacuum techniques for pressures ranging from 
300 mm down. Despite low luminosity, it was 
established that the filamentary spark transition 
from corona to a transient arc, as in air, pro- 
gressed from anode to cathode by means of pri- 
mary and secondary streamers down to 50 mm 

at which point the spark channel was diffuse with 
no recognizable streamers. Addition of 1% air 

at 60 mm of A restored the streamer sequence. 
Argon purified by prolonged gas cataphoresis 
showed no streamer mechanism at 240-mm pres- 
sure. Thus in the absence of adequate photoioni- 
zation by its own photons, the transition in pure 
A consists in the development of a diffuse arc- 
like channel on a time scale of microseconds in 
contrast to tenths of a microsecond with streamers 
when photoionizable impurities are present. 


*Now at Kyungpook National University, Taegu, Korea. 

tNow at Philips Laboratories, Irvington-on-Hudson, 
New York. 

tNow at University of Malaya, Kuala Lumpur, Malaya. 


ANISOTROPIC FERMION GAS. G. Wentzel, 
Enrico Fermi Institute for Nuclear Studies, 
University of Chicago, Chicago, Illinois (Re- 
ceived June 27, 1960). 


An energy term discovered by Kohn and Lut- 
tinger is shown to have a simple geometrical 
meaning. 


THERMAL CONDUCTION IN ROTATING LIQUID 
HELIUM II. C. E. Chase, Lincoln Laboratory, 
Massachusetts Institute of Technology, Lexing- 
ton, Massachusetts (Received June 30, 1960). 


When a heat current flows in a wide channel 
filled with liquid helium II, the resulting temper- 
ature gradient is approximately proportional to 
the cube of the heat current density. The estab- 
lishment of this gradient requires a time 7 which 
isa function of heat current, temperature, and 


past history of the helium. The present experi- 
ments concern the effect of uniform slow rota- 
tion, about an axis normal to the direction of 
heat flow, upon gradT and 7. GradT was meas- 
ured at 0 and 1.3 rad/sec, the highest angular 
velocity at which this measurement could con- — 
veniently be carried out. No effect of rotation 
could be observed; however, an approximate 
calculation suggests that grad7T might increase 
detectably at somewhat higher rates of rotation. 

T was measured at a number of angular velocities 
between 0 and 4 rad/sec; it was found that 7 was 
appreciably reduced by rotation, the effect being 
greatest at small heat currents and high angular 
velocities. These results can be explained on 

the assumption that mutual friction results from 
turbulence in the superfluid component, taking 
the form of a tangled mass of vortex line. The 
delay time 7 characterizes the rate of growth 

of this turbulence when a heat current is switched 
on; rotation reduces 7 by introducing an initial 
length of vortex line which accelerates this growth. 


IRREVERSIBLE PROCESSES IN ISOLATED SYS- 
TEMS. Edward N. Adams, IBM Research Center, 
Yorktown Heights, New York (Received May 10, 
1960). 


The point of view is taken that the irreversibility 
paradox of Loschmidt and Zermelo arises because 
there are two valid transport equations, one caus- 
al, the other anticausal, each consistent with the 
fundamental equations of mechanics. From this 
point of view the problem of irreversibility is to 
characterize a given nonequilibrium distribution 
as to which transport equation it will obey. Using 
the transport theories of Kohn and Luttinger and 
of Van Hove, we obtain a statistical criterion 
capable of so characterizing distribution functions 
of both stationary and time-varying types. We 
discuss the question of how experimental proce- 
dures consistently bring about a situation to which 
the causal transport equation applies; we answer 
the question for a very simple kind of experi- 
mental situation. 


SPIN SUSCEPTIBILITY OF AN ELECTRON GAS. 
P. A. Wolff, Bell Telephone Laboratories, Mur- 
ray Hill, New Jersey (Received June 24, 1960). 


The generalized random phase approximation 
is used to investigate the effects of electron- 
electron interaction on the spin susceptibility 
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of an electron gas. Within this approximation, 
the induced spin density is determined by the 
solution of an integral equation that describes 
exchange scattering between virtually excited 
particles. Solutions to this equation are obtained 
in simple cases. They indicate that the exchange 
scattering enhances the susceptibility over that 
predicted by the Ruderman-Kittel formula by 
factors which, for electron densities comparable 
to those found in good metals, are in the range 
1-2. These results can also be derived by sum- 
ming certain classes of diagrams in many-body 
perturbation theory. The structure of these 
graphs shows that the random phase approxima- 
tion takes account of the simplest sort of self- 
energy effect and the simplest type of scattering 
of virtually excited particles. 


ELECTRO-OPTIC KERR EFFECT AND POLAR- 
IZATION REVERSAL IN DEUTERIUM-DOPED 
ROCHELLE SALT. H. H. Wieder and D. A. 
Collins, U. S. Naval Ordnance Laboratory, 
Corona, California (Received June 30, 1960). 


Polarization reversal as a function of nuclea- 
tion and growth of domains in ferroelectric 
deuterium-doped Rochelle salt was investigated 
by means of the electro-optic Kerr effect. The 
results indicate that a phenomenological model 
based on statistical nucleation of domains in a 
plane including the ferroelectric axis followed 
by a two-dimensional sidewise expansion of the 
domains, will adequately account for the experi- 
mental observations for fields larger than 50 
volts/em. For lower fields, the process is con- 
trolled primarily by the nucleation of new do- 
mains due to localized stresses which hinder 
the displacement of domain walls. 


VACANCY-VACANCY INTERACTION IN COP- 
PER. V. G. Weizer and L. A. Girifalco, Lewis 
Research Center, National Aeronautics and 
Space Administration, Cleveland, Ohio (Received 
June 16, 1960). 


The binding energy of two vacancies in a static 
lattice as a function of their separation and the 
positions of their displaced neighboring atoms 
has been calculated using a Morse potentiai 
function model for copper. It was found that two 
vacancies attract one another at separation less 
than about 7 A. At separations greater than 7A 
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the vacancies do not interact appreciably. The 
most stable separation was found to be the first 
nearest neighbor separation or the divacancy 
configuration, for which the binding energy was 
found to be 0.64 ev. Based on these calculations, 
it is shown that third-stage annealing in irradiated 
copper may be accounted for by divacancy migra- 
tion. The role of the divacancy in copper self- 
diffusion is also explained. 


DE HAAS-VAN ALPHEN EFFECT IN GRAPHITE 
BETWEEN 3 AND 85 KILOGAUSS. W. J. Spry 
and P. M. Scherer, Research Laboratory, 
National Carbon Company, Parma, Ohio (Re- 
ceived June 22, 1960). 


The de Haas-van Alphen effect has been meas- 
ured in a single crystal of graphite between 3 
and 85 kilogauss with the c axis of the crystal 
parallel to the magnetic field. Two periods were 
observed: 2.05x10~° gauss™ for the electronic 
component and 1.53 x107° gauss~! for the holes. 
No new and unexpected oscillations were ob- 
served at magnetic fields as large as 100 kilo- 


gauss. 


FLUORESCENCE AND REFLECTION SPECTRA 
OF NaI SINGLE CRYSTALS. W. J. Van Sciver, 
Levinthal Electronic Products, Incorporated, 
Palo Alto, California (Received April 18, 1960). 


The fluorescence emission and excitation, re- 
flection, and transmission spectra of 22 different 
Nal crystal samples, grown in crucibles of Pt or 
quartz under atmospheres of N, or H, or H, plus 
L, and in some cases containing stoichiometric 
excess of Na or L, or containing various concen- 
trations of Tl, have been explored from 2.5 to 8 
ev. At -190°C an emission band centered at 295 
my (4.20 ev) was observed in all but the one 
sample which contained the highest Tl concen- 
tration, 0.1 mole %. Another emission band at 
375 mu (3.31 ev) was found to be strongest in 
crystals containing excess I,. A 425-my (2.92-ev) 
emission appeared strongly in all crystals con- 
taining Tl, but a 325-my (3.82-ev) band appeared 
only in the one crystal with 0.1 mole % Tl. With 
excitation energy below 5.5 ev the different emis- 
sion bands arise from discrete excitation bands. 
With excitation energies above 5.5 ev up to 8 ev 
the emission was primarily in the 295-my band 
until Tl concentration reached 0.05 mole % where 
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the 425-my band and 295-my band became ap- 
proximately equal. Maxima in reflection spectra 
coincide closely, but not exactly, with published 
maxima in the absorption spectra of NaI evapo- 
rated films. An inflection, not seen in the pub- 
lished absorption spectra, appears at 5.62 ev in 
the reflection spectrum at -190°C. Coincidental 
with the lowest energy maximum in the reflection 
spectrum at 5.59 ev, the excitation spectra for 
poth 295-my and 425-my emission exhibit minima 
too strong to be accounted for by reflection loss. 
At excitation energies of 5.7 ev and above, lumines- 
cence quantum efficiency is greatest for an illumi- 
nating angle of incidence of 45°, whereas at lower 
energies no definitely preferable incidence angle 
was found. Also, with higher exciting energies 

it was found that the emission was strongest on 

the axis of illumination. With exciting energies 
below 5.5 ev the emission was isotropic. 


MAGNETIC SUSCEPTIBILITY OF TETRAGONAL 
TITANIUM DIOXIDE. Frank E. Senftle, U. S. 
Geological Survey, Washington, D. C., Titus 
Pankey, Howard University, Washington, D. C. 
and Frank A. Grant, National Bureau of Standards, 
Washington, D. C. (Received July 1, 1960; re- 
vised manuscript received April 21, 1960). 


Careful measurements have been made of the 
magnetic susceptibility of the rutile and anatase 
crystalline forms of titanium dioxide. The mag- 
netic susceptibility of a single crystal of high- 
purity rutile was found to be (0.067+ 0.0015) x10~° 
emu per gram, and was temperature-independent 
from 55° to 372°K. Difficulty was encountered in 
obtaining a good value of the magnetic suscepti- 
bility of anatase because of impurities. However, 
a value of 0.02 x10~° emu per gram was obtained 
as a maximum value. A discussion is given for 
the different values obtained for anatase and 
rutile. 


PHOSPHORESCENCE MECHANISMS. I. AP- 
PROACH AND GENERAL ANALYSIS. C. Billing- 
ton, Division of Chemical Physics, Common- 
wealth Scientific and Industrial Research Organi- 
zation, Chemical Research Laboratories, Mel- 
bourne, Australia (Received June 13, 1960). 


A new technique is described for the study of 
decay mechanisms in phosphors. When a phos- 
phor is excited by radiation whose intensity is 


modulated to a depth less than 10%, the modula- 
tion of the emission intensity has negligible 
harmonic content. Quantitative comparisons of 
these two modulated signals, if of sufficient ac- 
curacy and range, result in qualitative distinc-. 
tions between different basic mechanisms, and 
the possibility of subsequent quantitative analy- 
sis. In particular, monomolecular and bimolec- 
ular processes are unequivocally distinguished. 


PHOSPHORESCENCE MECHANISMS. II. DE- 
SCRIPTION OF PHOSPHOROMETER. C. Billing- 
ton, Division of Chemical Physics, Common- 
wealth Scientific and Industrial Research Organi- 
zation, Chemical Research Laboratories, Mel- 
bourne, Australia (Received June 13, 1960). 


The functions of the phosphorometer are (i) to 
excite the phosphor under test with ultraviolet 
radiation whose intensity is modulated sinusoid- 
ally, (ii) to measure the modulation frequency, 
and (iii) to sample the excitation and emission 
intensities, and to extract from their modulated 
components the in-phase attenuation factor be- 
tween them. The Hg arc source is current- 
modulated and stabilized by optical feedback 
over the range 3 to 30000 cps. The average 
intensity of both excitation and emission is used 
to control the sensitivity of their respective 
photodetectors so that the modulation signals 
are always referred to the same datum. Both 
signals are passed in turn through a homodyne 
detector whose phase is adjusted to identity with 
the excitation modulation. The in-phase attenua- 
tion factor which is the ratio of the homodyne 
outputs is determined potentiometrically with a 
probable error of about 0.002. 


PHOSPHORESCENCE MECHANISMS. III. 
METHOD OF ANALYSIS. C. Billington, Division 
of Chemical Physics, Commonwealth Scientific 
and Industrial Research Organization, Chemical 
Research Laboratories, Melbourne, Australia 
(Received June 13, 1960). 


The method of deconvolution is discussed for 
the case in which phosphorescence is distributed 
between a number of concurrent components with 
different rate constants. The resolution of the 
method is deduced quantitatively from the known 
rms difference of the apparatus. 
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PHOSPHORESCENCE MECHANISMS. IV. DE- 
CAY RATE SPECTRA OF RUBY, URANIUM 
GLASS, AND Mn-ACTIVATED ZINC SULFIDE. 
C. Billington, Division of Chemical Physics, 
Commonwealth Scientific and Industrial Research 
Organization, Chemical Research Laboratories, 
Melbourne, Australia (Received June 13, 1960). 


Ruby, uranium glass, and a ZnS:Mn sample 
were studied by means of the phosphorometer 
in a program of measurement and analysis which 
included variation of temperature and average 
exciting intensity. It was not possible to identify 
any laws of temperature dependence, but there 
is strong evidence that part of the phosphores- 
cence mechanisms operating in the ZnS:Mn 
sample involves recombination. 


THEORY OF IMPURITY RESISTANCE IN METALS. 


J. S. Langer, Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania (Received June 27, 
1960). 


A many-body technique is developed for the 
calculation of the de resistivity of a Fermi fluid 
in the presence of a few, randomly scattered, 
fixed, impurities. A certain class of graphs 
yields an expression for the conductivity which 
is similar in form to the standard classical 
transport coefficient; but the decay time is de- 
termined by the scattering of single-particle- 
like excitations at the Fermi surface by screened 
impurities. A propagator method similar to that 
used in field theory is employed throughout the 
paper, and the perturbation-theoretic interpreta- 
tion of this method is examined in some detail. 


IMPURITY CONDUCTION AT LOW CONCEN- 
TRATIONS. Allen Miller* and Elihu Abrahams, 
Physics Department, Rutgers University, New 
Brunswick, New Jersey (Received June 23, 
1960). 


The conductivity of an n-type semiconductor 
has been calculated in the region of low temper- 
ature T and low impurity concentration Np: The 
model is that of phonon-induced electron hopping 
from donor site to donor site where a fraction 
K of the sites is vacant due to compensation. 

To first order in the electric field, the solution 
to the steady state and current equations is 
shown to be equivalent to the solution of a linear 
resistance network. The network resistance is 


220 


evaluated and the result shows that the T depend- 
ence of the resistivity is p <exp(¢,/kT). For 
small K, €, =(e?/k,)(4mmp/3)*(1 -1.35K™), where 
K, is the dielectric constant. At higher K, «, 
and p attain a minimum near K=0.5. The de- 
pendence on 7p is extracted; the agreement of 







the latter and of €, with experiment is satisfactory, | 


The magnitude of p is in fair agreement with ex- 
periment. The influence of excited donor states 
on p is discussed. 


*Present address: Physics Department, University 
of Illinois, Urbana, Illinois. 


FRICTIONAL STRESS ACTING ON A MOVING 
DISLOCATION IN AN OTHERWISE PERFECT 
CRYSTAL. Doris Kuhlmann-Wilsdorf, School 
of Metallurgical Engineering, University of 
Pennsylvania, Philadelphia, Pennsylvania (Re- 
ceived June 27, 1960). 


The problem of the frictional stress suffered 
by moving dislocations in otherwise perfect 
crystals is investigated. This is done without 
calculating the core energies of dislocations, 
but by considering stresses and strains on the 
slip plane. The level of frictional stresses ob- 
tained is much higher than reported previously. 
Since common glide dislocations in metals with 
close-packed structures apparently do not suffer 
significant frictional stresses, mechanisms are 
discussed which tend to reduce their effect. A 
new such mechanism is discovered. It is based 
on the idea that the positions of dislocation axes 
are not defined with precision, but only within 
one to a few times the average displacement of 
the oscillating atoms. The expected result of 
this, is a depression of the frictional stress for 
close-packed metals even at very low tempera- 
tures, almost no effect on dislocations in crys- 
tals with diamond structures, and a temperature 
dependence proportional to e~°°MStT/TM for 
NaCl type salts and, probably, for bcc metals. 


RELATION BETWEEN DIRAC AND CANONICAL 
DENSITY MATRICES, WITH APPLICATIONS TO 
IMPERFECTIONS IN METALS. N. H. March 

and A. M. Murray, Department of Physics, The 












University, Sheffield, England (Received June 
17, 1960). 


It is shown that the canonical density matrix 
in a single-particle framework may be related 












I 
t 
¢ 
i 
te 
8 
a 
nn 
of 
hi 
of 
th 
of 
inf 
loy 









sre 


tory. 


S 


e 
b- 
sly. 
ith 
iffer 


sed 
ixes 
‘in 

- of 
of 

3 for 
ra- 
rys- 
ature 
yr 
als. 


TICAL 
NS TO 
ch 

The 





‘une 


trix 
ated 





VoLUME 5, NUMBER 5 


PHYSICAL REVIEW LETTERS 





SEPTEMBER 1, 1960 





directly to the generalized canonical density 
matrix, containing the Fermi- Dirac function, 
and to the Dirac density matrix. 

A study is then made of density matrices in 
central field problems. A new differential equa- 
tion is derived, from the Bloch equation, for the 
diagonal element of the canonical density matrix. 
In the case of a continuum of energy levels, this 
is shown to lead directly to a differential equa- 
tion for the diagonal element of the Dirac ma- 
trix, that is the particle density. Free-electron 
density matrices are fully worked out and a per- 
turbation theory based on these free-electron 
forms is presented. 

It is further shown that for a nonspherical po- 
tential energy V(r), the work of Green on the 
quantum- mechanical partition function may be 
utilized to yield a perturbation theory for the 
Dirac matrix. In this way, the correct formula- 
tion to replace Mott’s well-known first order 
approximation for dealing with imperfections in 
metals is obtained. A brief discussion of the 
way in which this removes qualitatively the dif- 
ficulties of the Mott treatment is given, and the 
possibility of direct numerical application in a 
self-consistent framework is pointed out. 


DEPENDENCE OF SECONDARY ELECTRON 
EMISSION FROM MgO SINGLE CRYSTALS ON 
ANGLE OF INCIDENCE. A. B. Laponsky and 

N. R. Whetten, General Electric Research Labo- 
ratory, Schenectady, New York (Received 
November 24, 1959); revised manuscript received 
July 21, 1960). 


The influence of the angle of incidence of pri- 
mary electrons on the energy distribution and 
total yield of secondary electrons from MgO 
cleaved single crystals has been studied. Charg- 
ing effects were minimized by using single-pulse 
techniques. The relative number of low-energy 
secondary electrons decreases with increasing 
angles of incidence. The effect appears to be 
most marked at primary energies in the vicinity 
of the maximum of the yield curve. At lower and 
higher primary energies the dependence on angle 
of incidence diminishes. Under most conditions 
the total yield increases with increasing angles 
of incidence. The magnitude of the increase is 
influenced by the primary electron energy. At 
low energies a slight decrease in total yield with 
increasing incidence angles has been observed. 

A dependence of the backscattered fraction of 


electrons on the angle of incidence of the pri- 
maries has been observed at several primary 
electron energies. 

A number of maxima and minima occur super- 
posed on the gross variation of total yield with 
angle of incidence. The angular positions of the 
maxima are insensitive to primary energy but 
show a dependence on the crystal orientation. 
Several models for accounting for the phenomena 
are discussed. 


MECHANISM FOR COLORATION OF ALKALI 
HALIDES AT LOW TEMPERATURES. Clifford 
C. Klick, U. S. Naval Research Laboratory, 
Washington, D. C. (Received June 23, 1960). 


A mechanism for low-temperature coloration 
of alkali halides is proposed in which double 
ionization of a halide ion or single ionization of 
two adjacent halide ions leads to the formation 
of an uncharged halogen molecule which occupies 
two adjacent halide ion sites. If the molecule is 
localized at one of these sites by the jump of a 
nearby halide ion into the other site, a vacancy 
and interstitial are formed. Capture of an elec- 
tron by each of these imperfections leads to the 
production of the F center and H center. 


NUCLEAR QUADRUPOLE RESONANCE IN SU- 

PERCONDUCTING GALLIUM. R. H. Hammond 

and W. D. Knight, Department of Physics, Uni- 

versity of California, Berkeley, California (Re- 
ceived June 13, 1960). 


The nuclear quadrupole resonance (NQR) of 
Ga® is investigated between 0.8 and 4.2°K. (The 
superconducting critical temperature is 1.084°K.) 
Progressive saturation of the resonance is pro- 
duced at increasing rf power levels of a fre- 
quency- modulated marginal oscillator. The nu- 
clear spin-lattice relaxation rate 1/T, (sec™*) is 
approximately $7 (°K) in the normal state. The 
contact part of the hyperfine interaction appears 
to be predominant in producing relaxation. Al- 
though the saturation method does not permit a 
precise determination of the constant of propor- 
tionality, the comparison of normal and super- 
conducting relaxation rates is considerably more 
reliable. The rf fields are assumed to obey the 
London law of penetration, and an average rf 
field is determined in the superconductor. The 
relaxation rate shows a maximum enhancement 
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by a factor 1.8 at approximately 0.957,, a re- 
sult which agrees with that obtained by Slichter 
and Hebel, and Redfield for the nuclear magnetic 
resonance of aluminum, and which serves as an 
additional experimental justification for certain 
features of the Bardeen-Cooper-Schrieffer 
theory of superconductivity. Unsaturated signal 
intensities in normal and superconducting states 
furnish a basis for estimating the penetration 
depth in superconducting metal spheres: ) (0) 

= 1200 A for an average particle diameter of 

2.7 u. The NQR frequency of Ga in the super- 
conductor shifts by (+)510~° of the normal fre- 
quency [the corresponding result for indium is 
found to be approximately (-)10~?]. This means 
that the contribution of the conduction electrons 
to the average nuclear quadrupole coupling is 
modified by the rearrangement of the conduction 
band in the superconductor. It is demonstrated 
that, if the quadrupole term in the hyperfine in- 
teraction were to predominate, the spin relaxa- 
tion rate in the superconductor would have a 
temperature dependence like that of the ultra- 
sonic attenuation. 


NEW METHOD FOR TREATING THE ANTI- 
FERROMAGNETIC GROUND STATE. H. L. 
Davis,* Department of Physics and Astronomy, 
The Ohio State University, Columbus, Ohio (Re- 
ceived January 25, 1960; revised manuscript 
received June 22, 1960). 


A perturbation expansion for the ground-state 
energy of an antiferromagnetic spin system is 
obtained in terms of a linked-spin-cluster ex- 
pansion similar to Goldstone’s linked- Feyman- 
diagram expansion for the interacting Fermion 
system. From the energy perturbation series, 
perturbation series for the long- and short-range 
order may be obtained. Using these perturbation 
series, the ground state properties are calcu- 
lated through seventh order and compared with 
the results obtained by other investigators. In 
all cases, the values obtained here for the 
ground state energy are lower than those which 
have been obtained by purely variational means. 
The results for the long-range order are radi- 
cally different from the variational results but 
agree qualitatively with those obtained by spin- 
wave theory; however, the method is free of the 
usual objections which are voiced to spin-wave 
treatments of antiferromagnetism. The present 
work is incomplete in that limits on the error 
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introduced by using only a finite number of terms 
of the perturbation series to calculate the physi- 
cal properties are not obtained. But the author 
feels that the merit of the present work is in the 
method rather than the results since it provides 
a consistent framework both to settle the con- 
vergence question and to treat the antiferromag- 
netic spin system at finite temperatures. 


*Now at Physical Sciences Research Department, 
Sandia Corporation, Albuquerque, New Mexico. 


VALLEY-ORBIT SPLITTING OF ANTIMONY 
IN GERMANIUM. H. Fritzsche, Department of 
Physics and Institute for the Study of Metals, 
University of Chicago, Chicago, Illinois (Re- 
ceived June 23, 1960). 


The change of electrical conductivity under 
uniaxial tension and compression has been meas- 
ured over the range 4°K to 7°K for single-crystal 
specimens of germanium doped with antimony. 
The stress was varied from 1 x10" to 5x10 
dynes/cm?. On the basis of Price’s calculation 
of the effect of shear on the Kohn-Luttinger donor 
level structure, an expression for the piezoresis- 
tance has been derived, which includes terms 
of high order in the strain. It is shown that for 
a finite valley-orbit splitting, i.e., a finite energy 
separation between the onefold and the threefold 
1s-like donor states, shear increases the total 
electron concentration in the conduction band. 
For uniaxial stresses along the [110] direction 
this increase in electron concentration is an even 
function of stress and can, therefore, be de- 
termined from a linear combination of the piezo- 
resistance measured under tension and com- 
pression. A comparison of the theoretical ex- 
pression with the experimentally obtained change 
in electron concentration yields for antimony in 
germanium a valley-orbit splitting of 0.57 + 0.03 
milli-electron volt. 


QUANTUM MECHANICAL LAW OF CORRE- 
SPONDING STATES FOR VAN DER WAALS 
SOLIDS AT 0°K. Newton Bernardes, Institute 
for Atomic Research and Department of Physics, 
Iowa State University, Ames, Iowa (Received 
June 22, 1960). 


A quantum mechanical variational treatment 
based on a simple Heitler-London wave function 
is used to describe the various properties of 
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Van der Waals solids at 0°K. The effects of nu- 
clear motion on the cohesive energy, volume, 
compressibility, sound velocity, etc., are dis- 
cussed. These effects can be expressed in the 
form of a power series inf, of which we show 
the first term, linear in #, explicitly. The re- 
sults are in good agreement with the available 
experimental data for all solidified inert gases, 
except for He to which the present method is not 
applicable. 


ELECTRICAL RESISTIVITY OF THE HEAVY 
RARE-EARTH METALS. R. V. Colvin, * Sam 
Legvold, and F. H. Spedding, Institute for Atomic 
Research and Department of Physics, Iowa State 
University, Ames, Iowa (Received June 8, 1960; 
revised manuscript received August 1, 1960). 


The electrical resistivities of polycrystalline 
Gd, Tb, Dy, Ho, Er, Tm, and Lu have been 
measured between 1.3°K and room temperature. 
The slope of the resistivity curve for Gd changes 
near the Curie point. The curve for Tb is very 
much like that for Gd but there is some evidence 
that two ordering temperatures exist for this 
metal. Dy, Ho, and Tm all show peaks in resis- 
tivity near their Néel points, while Er shows only 
a change in slope at its Néel point. The change 
from ferromagnetism to antiferromagnetism in 
Dy is seen as a sharp rise in the resistivity. 


*Now at the Edgar C, Bain Laboratory of the U. S. 
Steel Corporation, Monroeville, Pennsylvania. 


EFFECT OF POINT IMPERFECTIONS ON LAT- 
TICE THERMAL CONDUCTIVITY. Joseph Call- 
away, Division of Physical Sciences, University 
of California, Riverside, California, and Hans C. 
von Baeyer, Department of Physics, University 
of Miami, Coral Gables, Florida (Received June 
1, 1960). 


The consequences of a simple, phenomenological, 
theory of lattice thermal conductivity with respect 
to the effect of point imperfections are summar- 
ized. The experimental results of Berman et al. 
on the effect of varying the concentration of Li® 
on the conductivity of lithium fluoride are ana- 
lyzed in detail. 


THEORY OF THE COMPLEX REFRACTIVE 
INDEX. C. Alden Mead, School of Chemistry, 
University of Minnesota, Minneapolis, Minnesota 
(Received February 12, 1960). 


A quantum-mechanical formalism is developed 
for representing the complex refractive index of 
a gas, including both dispersive and absorptive 
parts. In analogy with some previous work by 
the author, it is shown that the phenomenological 
Maxwell equations are obeyed as equations be- 
tween certain matrix elements of operators rep- 
resenting the field quantities. The method makes 
use of a canonical transformation due to Arnous 
and Bleuler, and also of the “damping operators” 
previously introduced by the author. It is shown 
that the calculation of the complex refractive 
index reduces to the calculation of a single 
diagonal matrix element of a damping operator. 
The actual evaluation of this matrix element is 
taken up in the following paper. 


RESONANCE BROADENING OF SPECTRAL 
LINES. C. Alden Mead, School of Chemistry, 
University of Minnesota, Minneapolis, Minnesota 
(Received February 12, 1960). 


A formalism introduced previously by the 
author is used to calculate the shape of an ab- 
sorption line in a gas, taking account of the 
“resonance” interaction between like atoms. 
The most serious approximation is the neglect 
of the effect of translational motion. An asym- 
metric line is obtained, with the maximum 
shifted slightly toward the red. Some simple 
classical considerations also agree qualitatively 
with the results. Various properties of the line 
shape obtained are briefly discussed. The neg- 
lect of translational motion makes it impossible 
to compare the results with experiments per- 
formed up to the present time. 


ELASTIC AND INELASTIC SCATTERING OF 
ELECTRONS FROM THE S STATES OF ATOMIC 
HYDROGEN. Kenneth Smith, Argonne National 
Laboratory, Argonne, Illinois (Received May 19, 
1960; revised manuscript received August 1, 
1960). 


Total elastic and inelastic cross sections, for 
the scattering of electrons (including exchange) 
from the ls, 2s, and 3s states of atomic hydro- 
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gen, have been calculated in three approxima- 
tions using the method of partial waves. The 
purpose of this report is to present the S-, P-, 
and D-wave contributions to the elastic, Q1.>_15, 
and inelastic, Q;,_9,, cross sections. 


METASTABLE HYDROGEN MOLECULES. 
William Lichten, Ryerson Laboratory, Depart- 
ment of Physics, University of Chicago, Chicago, 
Illinois (Received June 30, 1960). 


The existence of a metastable electronic state 
of H, has been established by means of the molec- 
ular beam magnetic resonance technique. Ob- 
servations have been made of the radio-frequency 
Zeeman spectrum of natural and para-hydrogen. 
At a magnetic field of about 5 gauss the strongest 
line corresponds to a g value of 0.83. The ap- 
pearance potential for this line is 11.9+0.25 ev. 
Observation of the quadratic Zeeman effect made 
it possible to estimate the fine structure of the 
N =2 level of para-hydrogen. The observed in- 
tervals are: (J=2—J=3)=5730+ 200 Mc/sec, 

(J =2-— J=1)=5000+ 200 Mc/sec. This repre- 
sents the first complete resolution of the fine 
structure of a rotational level of H,. Knowledge 
of the electronic structure of H,, agreement of 
the appearance potential with spectroscopic data, 
agreement of the fine structure with partially 
resolved optical spectra, and good agreement of 
the Zeeman spectra with theory lead to the con- 
clusion that the metastable state must be H,, 
c(log, 17,,) *T,. 


HYPERFINE STRUCTURE OF HYDROGEN, DEU- 
TERIUM, AND TRITIUM. L. Wilmer Anderson* 
and Francis M. Pipkin, Lyman Laboratory, 
Harvard University, Cambridge, Massachusetts, 
and James C. Baird, Jr.,f Mallinckrodt Labo- 
ratory, Harvard University, Cambridge, Mass- 
achusetts (Received June 3, 1960). 


The optical transmission of an optically oriented 
rubidium vapor in spin-exchange equilibrium with 
atomic hydrogen, deuterium, and tritium has been 
used to measure with high precision the hyperfine 
splittings of these paramagnetic atoms. The re- 
sults are Av(H) = 1420.405726 + 0.000030 Mc/sec, 
Av(D) = 327.384349 + 0.000005 Mc/sec, and d(T) 
= 1516.701396+ 0.000030 Mc/sec. These results 
are based on a value of the hyperfine splitting of 
Cs*** which is taken to be Av(Cs**5) = 9192.631840 
Mc/sec. 
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These measurements were made in various 
buffer gases which caused a shift in the observed 
hyperfine splitting, and the results given repre- 
sent extrapolations to zero pressure. The pres- 
sure shifts were measured for H in argon, neon, 
helium, and molecular hydrogen and were meas- 
ured for D and T in argon and neon. The assigned 
limits of error represent the range of disagree- 
ment of the zero-pressure extrapolations in the 
different buffer gases. 


*Now in the Physics Department, University of 
Wisconsin, Madison, Wisconsin. 

t+Now at the California Research Corporation, Rich- 
mond, California. 


APPLICABILITY OF ROOTHAAN’S SELF-CON- 
SISTENT FIELD THEORY. S. Huzinaga,* Labo- 
ratory of Molecular Structure and Spectra, De- 
partment of Physics, University of Chicago, 
Chicago, Illinois (Received June 28, 1960). 


Roothaan’s open-shell self-consistent field 
theory is extended in such a way that it can cover 
a number of important classes of atomic and 
molecular states which are not included in the 
range of applicability originally claimed. 


*Permanent address: Kyushu University, Fukuoka, 
Japan. 


POSITRONIUM DECAY IN MOLECULAR SUB- 
STANCES. Werner Brandt, Radiation Physics 
Laboratory, Engineering Department, E. I. 

du Pont de Nemours and Company, Incorporated, 
Wilmington, Delaware, and S. Berko and W. W. 
Walker, Department of Physics, University of 
Virginia, Charlottesville, Virginia (Received 
March 16, 1960). 


The overlap of the positron component of the 
Swiss-cheese wave function of positronium (Ps) 
in dense molecular substances with the lattice 
wave function determines the 2y annihilation 
rate of positrons bound in orthopositronium via 
electron pickup from the lattice (7, decay), and 
in part the 2y angular correlation. Ps acts as 4 
probe for the dependence of one-particle wave 
functions in lattices on structural lattice para- 
meters, and of radical-molecule interactions on 
atomic lattice parameters. To explore the sensi- 
tivity of this probe Ps wave functions have been 
calculated for different lattice structures in the 
Wigner-Seitz approximation, and the dependence 
of 7, decay on lattice parameters derived. This 
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dependence was tested experimentally by 7, meas- 
urements on carefully characterized molecular 
substances in different physical states. It is 
confirmed that the temperature and state de- 
pendence is primarily a “free volume” effect in 
the sense that the overlap integral between the 
Ps and lattice wave functions decreases with 
increasing lattice spacing. The effects of a 
spherical confinement on the self-annihilation 
rates and stability of Ps are shown to be small 
for radii relevant to molecular lattices. 

An anomalous behavior is found for the ice- 
water transition, where T, increases significantly 
despite the contraction of water between the melt- 
ing point and 4°C. 


K-ELECTRON CAPTURE BRANCH OF Sr*””, 

A. W. Sunyar and M. Goldhaber, Brookhaven 
National Laboratory, Upton, New York (Received 
June 29, 1960). 


The existence of a (0.65+ 0.25)% K-electron 
capture branch in the decay of 2.8-hr Sr®’” has 
been established by means of an “internal source” 
technique. The low logft=4.25 and small energy 
release of 115 kev make this K branch especially 
interesting because it establishes Rb®’ as a pos- 
sible low-threshold detector for the inverse 
neutrino capture reaction, Rb®’+v—Sr®’” +e-, 
which could be identified by observing the char - 
acteristic decay of the Sr®’”™ isomer. 


RADIOACTIVE DECAY OF Lu’. R. G. Wilson 
and M. L. Pool, Department of Physics and 
Astronomy, Ohio State University, Columbus, 
Ohio (Received April 4, 1960). 


Ytterbium oxide enriched to 81.4% in the 170 
mass number was irradiated with 6- Mev protons. 
An activity decaying by electron capture with a 
half-Jife of (2.054 0.05) days was produced and 
assigned to Lu’”° by the identification of the ytter- 
bium K x ray and by comparison with the activi- 
vities produced by similar proton irradiations of 
each of the other enriched isotopes of ytterbium. 
The observed activity consists of the ytterbium 


‘Kx ray, gamma rays with energies of 84, 193, 


245, 1010, 1080, 1165, 1275, 1415, 2035, 2365, 
2665, 2890, and 3085 kev, and a small amount of 
annihilation radiation. Because no particle radia- 
tion exists, the mode of decay of Lu’” is solely 
by electron capture to Yb'”°. Gamma-gamma 


coincidence measurements have led to the postu- 
lation of levels in Yb*”° at 2120 (0), 2365 (1 or 2), 
3170 (3), and 3395 (1 or 2) kev in addition to the 
previously assigned 84 (2+), and 278 (4+) levels. 
A partial energy level scheme with approximate 
electron-capture branching ratios is proposed. 


APPROXIMATE CALCULATION OF NUCLEAR 
BINDING ENERGY. Leonard S. Rodberg and 
Vigdor L. Teplitz, University of Maryland, 
College Park, Maryland (Received June 16, 
1960). 


The binding energy, size, and shape of finite 
nuclei are studied using a method which tests 
the “local uniformity” assumption. This assump- 
tion implies that the properties of finite nuclei 
can be obtained from those of infinite nuclear 
matter. The energy is computed using as a 
trial wave function an amplitude- and frequency- 
modulated plane wave; this permits the use of 
nuclear -matter results in a straightforward 
manner. The calculation gives a binding energy 
which is too low, a radius which is too small, 
and a nuclear surface which is too diffuse. 
Several possible corrections to this result are 
examined. 


ATOMIC BEAM STUDY OF THE HYPERFINE 
STRUCTURE OF THULIUM- 170. Amado Y. 
Cabezas and Ingvar Lindgren, * Lawrence Radi- 
ation Laboratory, University of California, 
Berkeley, California (Received June 3, 1960). 


The atomic-beam magnetic-resonance technique 
has been used to measure atomic and nuclear 
quantities of the 129-day isotope Tm?” as follows: 
J=7/2,g7=1.14122+ 0.00015, J=1, |A| =20043 
Mc/sec, and |B| =1010+15 Mc/sec. The values 
of J and g, are consistent with the ground-state 
assignment *F,,.. Values of the nuclear moments 
are calculated from the hyperfine- structure inter- 
action constants A and B by use of a two-para- 
meter radial wave function, in which one para- 
meter is determined from comparison with 
Hartree functions and the other parameter from 
the experimental spin-orbit coupling constant. 
Uncorrected values are obtained, as follows: 
| yy! =0.26+ 0.02 nm and |Q| =0.61+ 0.05 barn, 
with the two moments of the same sign. The same 
wave function is used to calculate the relativistic 
and diamagnetic corrections to the atomic g value, 
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and the result is in excellent agreement with the 
experiment. 


*On leave from the Institute of Physics, University 
of Uppsala, Uppsala, Sweden. 


PREFORMATION FACTOR IN EMISSION OF 
COMPLEX PARTICLES FROM NUCLEAR RE- 
ACTIONS. Bernard L. Cohen, University of 
Pittsburgh, Pittsburgh, Pennsylvania (Received 
May 27, 1960; revised manuscript received July 
11, 1960). 


The resason why complex particles (alphas, 
tritions, Li® nuclei, etc.) are emitted in nuclear 
reactions as frequently as nucleons (after cor- 
rections for Coulomb barrier penetration and 
energetics), whereas ice crystals are never 
emitted from evaporating water droplets, is 
investigated. It is shown that the difference is 
entirely explained by the fact that a nucleus is a 
highly degenerate system subject to Fermi-Dirac 
statistics, whereas an evaporating water droplet 
is a nondegenerate statistical system. 


NEW NUCLIDE, KRYPTON-74. J. H. Gray, II, 
H. L. Hamester,* and A. A. Caretto, Jr., De- 
partment of Chemistry, Carnegie Institute of 
Technology, Pittsburgh, Pennsylvania (Received 
June 13, 1960). 


A new nuclide of krypton, Kr™, has been pro- 
duced by high-energy proton interaction with 
strontium. It was found to have a half-life of 
(20+1) minutes and to decay by (3.10+ 0.1)-Mev 
positrons. No gamma rays in the energy range 
0.4 to 4.0 Mev were observed. Evidence for the 
production of Kr” was also observed. Calcula- 
tions indicate an upper limit for the half-life of 
(0.8+ 0.4) minute. 


*Present address: Department of Chemistry, Uni- 
versity of New Mexico, Albuquerque, New Mexico. 


ASYMMETRY IN THE ELASTIC SCATTERING 
OF POLARIZED PROTONS BY DEUTERIUM. 
Stewart M. Marcowitz, The Enrico Fermi Insti- 
tute for Nuclear Studies and Physics Department, 
The University of Chicago, Chicago, Dlinois 
(Received June 3, 1960). 


Asymmetries in the elastic scattering of polar- 
ized protons by deuterium have been measured 
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at 419 Mev at proton center-of-mass angles 
from 47.9° to 115°. The inferred polarizations 
are: -0.092+ 0.058 at 47.9°, -0.502+ 0.058 at 
59°, -0.583+ 0.081 at 68°, -0.128+ 0.046 at 77.8°, 
-0.122+ 0.111 at 90°, -0.141+0.056 at 97°, -0.309 
+ 0.073 at 105°, and +0.174+ 0.059 at 115°. It is 
suggested that using the impulse approximation 
these results may serve as a test of nucleon- 
nucleon scattering phase shifts or, with proton- 
proton scattering data, may be used to obtain 
Z=0 nucleon-nucleon phase shifts. 


RARE £2 TRANSITION IN C*’. R. E. Pixley, 

J. V. Kane, and D. H. Wilkinson,* Brookhaven 
National Laboratory, Upton, New York (Received 
June 22, 1960). 


By means of a coincidence method it is shown 
that the 5/2+ state of C’* at 3.86 Mev has an E2 
branch of relative strength (9.3 + 2.0) x107° to 
the $+ state at 3.09 Mev. From consideration of 
the likely strengths of the competing £1 transi- 
tion to the 3/2- state at 3.68 Mev and the M2 
transition to the ground state, it is deduced that 
this branching ratio implies an E2 transition 
speed of the order of a single-particle (proton) 
unit. This in turn demands the substantial partici- 
pation of at least one excited state of C’ in the 
parentage of at least one of the C’® states with 
respect to the 1d or 2s neutron. This eliminates 
in particular a jj-coupling description of C** in 
which the two states in question are 1d, and 
2s, neutron states and more generally a weak- 
coupling model in which the ground state of C” 
is the unique parent for the C states (with re- 
spect to 1d and 2s neutrons). 


*Permanent address: Clarendon Laboratory, Oxford, 
England. 


MEASUREMENT OF THE O'°(n, p)N*® CROSS 
SECTION AT 14.7 Mev. J. A. DeJuren and R. W. 
Stooksberry, Bettis Atomic Power Laboratory, 
Westinghouse Electric Corporation, Pittsburgh, 
Pennsylvania (Received June 22, 1960). 


The O'*(n, p)N*® cross section has been meas- 
ured for 14.7-Mev neutrons. A recoil proton 
telescope was used to measure the neutron flux. 
Oxygen was incorporated in a liquid scintillator 
which was irradiated in the flux. The absolute 
activity of the resulting N’* was measured by 47 
beta counting the solution with two photomulti- 
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pliers in coincidence. A value of 39.2+1.6 mb 
was obtained for the cross section. 


EXTENSION OF THE SHELL MODEL FOR 
HEAVY SPHERICAL NUCLEI. Michel Baranger, 
Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received June 17, 1960). 


The Bardeen-Bogoliubov-Belyaev treatment of 
the pairing correlations is applied to spherical 
nuclei with a general nuclear force. The inter- 
action between quasi-particles is treated by the 
method of linearized equations of motion. An 
advantage of this treatment is that the same 
equations describe single-particle excitations 
and collective excitations, so that the former 
are orthogonal to the latter and the total num- 
ber of states is correct. Another advantage is 
that the spurious states due to the fluctuations 
in the number of particles are automatically 
eliminated. The equations to be solved resem- 
ble those for a two-body shell model calculation. 
Simple estimates, based on delta-function or 
quadrupole forces, are made for the vibrational 
frequencies in various modes and transition 
matrix elements. It is concluded that the method 
is as powerful as other known methods for deal- 
ing with collective states by the shell model, 
and that the same order of magnitude for the 
effective nuclear force seems capable of fitting 
all the data. 


COMPARISON OF Po”!° ALPHA- PARTICLE EN- 
ERGY WITH THE Li’(p,”)Be’ REACTION THRES- 
HOLD ENERGY. C. P. Browne, J. A. Galey, J. R. 
Erskine, and K. L. Warsh, Department of Phys- 
ics, University of Notre Dame, Notre Dame; 
Indiana (Received June 17, 1960). 


Recent absolute measurements of the Po?!° 
alpha-particle energy disagree with the older 
value used as the standard for many nuclear re- 
action energy measurements. A new comparison 
with the Li’(p,m)Be” reaction threshold energy 
was made using the Notre Dame electrostatic 
accelerator and broad-range spectrograph. Four 
separate methods of comparison were used. In 
the first three the threshold was run and then 
protons or deuterons were scattered from appro- 
priate targets so that the scattered group was 
recorded on the spectrograph plate near the alpha 
group from a source placed at the target posi- 


tion. First, the spectrograph and, second, the 
beam analyzer were used to compare particle 
momenta. Third, with both fields held constant 
after the threshold was run with the molecular 
beam, deuterons were scattered, giving par- 
ticles of the same Bp as the alphas. In the 

fourth method several reaction energies that are 
precisely known in terms of the Li’(p,n)Be’ re- 
action threshold energy were measured in terms 
of the Po*’° alpha-particle energy. These were 
the Mg**(d, d’)Mg*** reaction to the first excited 
state of Mg™ and the N**(d, p)N"® reaction leading 
to three excited states of N’®. The four measure- 
ments agree and give 5.3086+ 0.003 Mev for Po?!° 
alpha-particle energy based on 1.8811 Mev for 
the Li’(p,)Be’ reaction threshold energy. 


jj-COUPLING MODEL IN ODD-ODD NUCLEI. 

M. H. Brennan* AND A. M. Bernstein, Palmer 
Physical Laboratory, Princeton University, 
Princeton, New Jersey (Received June 10, 1960). 


An analysis has been made, using the jj-cou- 
pling model, of the spins of 76 low-lying levels 
of odd-odd nuclei with 20< A< 120. Levels have 
been excluded from the analysis if any ambiguity 
exists in the assignment of a configuration. Ex- 
cluding particle-hole configurations, we find 
Nordheim’s “strong” rule obeyed in the 22 cases 
in which it is applicable. Nordheim’s “weak” 
rule is replaced by a rule predicting a ground 
state of the highest or lowest allowable spin. 
This revised rule is obeyed in 38 out of 41 cases. 
This competition between two levels of greatly 
differing spin results in the frequent occurrence 
of isomerism. For particle-hole. configurations 
the state of the highest spin minus one is the 
ground state in 6 out of 13 cases. This agree- 
ment with experiment is obtained using proton 
and neutron configurations which, in 66 of the 
76 cases, are found as ground states in the 
neighboring odd-even nuclei. The three revised 
coupling rules are predicted by the calculations 
of Schwartz, in which the residual proton-neutron 
interaction has a delta-function radial dependence, 
if the singlet-to-triplet strength is taken as 0.6, 
independent of mass number. This value is the 
same as that required to fit the free two-nucleon 
data. Exceptions to the empirical coupling rules 
and this theory will be discussed. 


*Now at Project Matterhorn, Princeton University, 
Princeton, New Jersey. 
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RECOIL PROPERTIES OF FISSION PRODUCTS. 
John M. Alexander and M. Frances Gazdik, Law- 
rence Radiation Laboratory, University of Cali- 
fornia, Berkeley, California (Received May 16, 
1960). 


The range and relative rates of energy loss in 
Al and Au have been measured radiochemically 
for five products from thermal-neutron- induced 
fission of U***. Range-velocity relationships for 
the median light product of fission and the median 
heavy product have been obtained from these 
measurements and other workers’ energy-loss 
data. The relation of range (R) to energy (E) or 
velocity (V) can be fitted to functions of the form 
R=kV -A or R=KE**. We have assumed that 
these functional forms can be applied to fission 
products of any mass. The constants K and A 
were determined from values of the range and 
kinetic energy for products of high yield. The 
values of these constants have been extrapolated 
to products of low yield. We have estimated 
kinetic energies, heretofore unmeasured, from 
the ranges of low-yield products. 

We have interpreted certain radiochemical 
observations in terms of the average component 
of the range perpendicular to the original velocity. 
The value of this component in Au has been esti- 
mated to be about one-fifth the total range. 


ISOMERS OF Am”™?. F. Asaro, I. Perlman, 

J. O. Rasmussen, and S. G. Thompson, Lawrence 
Radiation Laboratory and Department of Chem- 
istry, University of California, Berkeley, Cali- 
fornia (Received June 9, 1960). 


New measurements on the Am™? isomers 
show that the 16-hr activity is the ground state 
and the 152-yr activity the upper isomer, con- 
trary to earlier evidence. By electrostatic 
collection techniques some 16-hr americium 
activity was separated from an aged sample 
containing the 152-yr isomer. Careful study 
of the conversion electron spectrum of the sam- 
ple containing 152-yr Am™?” revealed eleven 
conversion lines of a 48.6-kev transition con- 
verted in americium, and the relative M-sub- 
shell conversion intensities are only consistent 
with an E4 assignment. This information, coupled 
with the recent atomic beam measurement of spin 
1 for the 16-hr activity, leads to a spin assign- 
ment of 5 for Am*#?””"_ 

The experimental information on these iso- 
mers is given a detailed analysis in terms of 
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the Bohr-Mottelson model. The most surprising 
conclusion of this analysis is that the spin 1 
ground state has a K-quantum number of Zero, 
and hence a negative spectroscopic quadrupole 
moment. Data on levels and decay of the low- 
spin isomers of Ho’ and Ta’®° are also con- 
sidered and lead to the conclusion that there is 
generally an energy displacement between odd 
and even spin members of K =0 rotational bands 
in odd-odd spheroidal nuclei. 





NUCLEAR ORIENTATION OF Mn™ AND Mn*”” 
Rudolf W. Bauer, Martin Deutsch, Gordon S. 
Mutchler, and Donald G. Simons,* Department 
of Physics and Laboratory for Nuclear Science, 
Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts (Received June 17, 1960). 


The spins of the nuclides Mn™ (290 days) and 
Mn*”""(21 min), incorporated in the lattice of 
cerium magnesium nitrate and nickel fluosilicate 
crystals, have been oriented at low temperatures. 
The measurement of the angular distribution of 
the gamma rays from Mn™ indicate the spin of 
Mn™ to be 3 or 2; in the latter case, the beta 
decay must be predominantly of the Fermi type 
with an upper limit of 10% of a possible Gamow- 
Teller admixture. From a simultaneous meas- 
urement of the angular distribution of the gamma 
rays from Mn™ and 5.7-day Mn™ the ratio of 
the nuclear g values of the two isotopes is found, 
which yields a magnetic moment for Mn™ of 
2.55+ 0.21 nm (for spin 3) or 2.16+ 0.26 nm (for 
spin 2). The circular polarization measurements 
of the gamma rays from Mn™ determine this 
moment to be positive. Similar angular distri- 
bution experiments on Mn**”, assuming spin 2, 
yield a magnetic moment of 1.04+0.16 nm if 
the beta decay is predominantly a Gamow-Teller 
transition, or 0.72+0.16 nm for a pure Fermi 
transition. These results indicate that the nu- 
clear g values of the 21-min and the 5.7-day 
states of Mn™ are about the same. 


*On leave from U. S. Naval Ordnance Laboratory, 
White Oak, Silver Spring, Maryland. | 








ASYMMETRY IN 143-Mev pn SCATTERING. 
Arthur F. Kuckes* and Richard Wilson, Cyclotron 
Laboratory, Harvard University, Cambridge, 
Massachusetts (Received June 16, 1960). 


i i sii i 


A beam of 145-Mev protons of polarization 
69+ 2% was scattered from liquid deuterium. 
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pm scattering events were identified by a coinci- 
dence method. The asymmetry in the scattering 
was measured. These data can be compared with 
polarization in free mp scattering at the appro- 
priate angles. 


*Now at Project Matterhorn, Princeton University, 
Princeton, New Jersey. 


ELASTIC SCATTERING OF 5-Bev 17> MESONS 
ON HYDROGEN. Richard G. Thomas, Jr., 
Lawrence Radiation Laboratory, University of 
California, Berkeley, California (Received 
June 9, 1960). 


The elastic scattering of 5-Bev negative pions 
on hydrogen has been observed in a large propane 
bubble chamber. We report below the results of 
this observation. It is shown that ~=(7+3)% of 


the events called elastic are probably background. 


The observed angular distribution is strongly 
peaked forward. This fact suggests that diffrac- 
tion scattering is the dominant process at this 
energy. The theoretical analysis of the distri- 
bution is in terms of the optical model. It is 
shown that the proton acts like a partially opaque 
sphere of radius 1.04107 cm +5%. The total 
elastic cross section is found to be 5.6+0.5 mb. 
From the extrapolated value of do(0)/dQ =29.8 
mb/sr + 10% in the center-of-mass system, a 
value of 29.1+ 2.9 mb was calculated for the total 
hydrogen cross section. The opacity of the 
sphere is thus 0.69 + 0.05. 


INTERACTION OF K* MESONS WITH EMULSION 
NUCLEI; K*-ENERGY RANGE 140 - 375 Mev. 

B. Sechi- Zorn and G. T. Zorn, Brookhaven Na- 
tional Laboratory, Upton, New York (Received 
June 6, 1960). 


Experimental data are presented on the inter- 
action in nuclear emulsions of Kt mesons in the 
energy range 140 - 375 Mev. 

From an optical- model calculation based on a 
Wood-Saxon potential (7,= 1.154 f, d=0.57 f), 
the real K-nuclear potential V was found to be 
positive (or repulsive), varying from 18.5+ 3.3 
Mev at the lowest energy, to 13.5+ 5.0 Mev at the 
highest energy. The value of the imaginary nu- 
clear potential W changed with increasing energy 
from -12.9+ 1.4 Mev to -17.642.5 Mev. 

The K-nucleon cross section was derived from 
the imaginary nuclear potential and, under the 


assumption of the independent-particle model of 
the nucleus, the total K-neutron non-charge- 
exchange cross section ox, + and the charge- 
exchange cross section o Kno were determined. 

It was found that oy,+ decreases with increasing 
energy from 12.1+1.6 to 8.1+1.5 mb, whereas 
Oxyo increases from 4.6+ 1.2 mb to 10.34 2.4 mb. 

A phase-shift analysis using the total cross 
sections 0x + ( the K-proton cross section), 
OKno» Go, b (the average K-nucleon cross section), 
and V gave evidence for P-wave contribution in 
isotopic spin state T=0. Changes in the K-neutron 
integrated cross sections with increasing K* 
energy seem to be accompanied by a change from 
a predominantly forward or near symmetric toa 
predominantly backward center-of-mass angular 
distribution for the charge-exchange K-neutron 
scattering. The opposite follows for the non- 
charge-exchange K -neutron differential cross 
section. 

Charged 7-meson production by K* mesons has 
been observed, and the K-nucleon cross section 
at ~280 Mev K* energy was estimated to be ~1/7 
mb. 


EXPERIMENTAL STUDY OF THE MAGNETIC 
STRUCTURE OF THE NEUTRON. J. I. Fried- 
man, H. W. Kendall, and P. A. M. Gram, High- 
Energy Physics Laboratory, Stanford University, 
Stanford, California (Received June 16, 1960). 


A measurement of the ratio of the magnetic 
form factor of the neutron to that of the proton 
has been carried out by comparing large- and 
small-angle elastic electron-deuteron scattering 
at constant four-momentum transfers. The ex- 
perimental result for the average value of the 
ratio in the range of momentum transfers from 
1.6 f= to 2.25 f™ is Fi,/F » = (0.91 + 0.05) + 0.07; 
the first error is a standard deviation from ex- 
perimental uncertainties, and the second from 
theoretical uncertainties in the analysis. Meas- 
urements of the ratio of the nucleon isotopic 
scalar form factors have also been obtained 
from this experiment. The average value of 
F,5/F,5 for the same range of momentum trans- 
fers has been found to be (+0.06+ 0.09) + 0.13. 
The small-angle scattering data have been used 
to determine the charge form factor of the deu- 
teron in the range of momentum transfers from 
0.98 f-" to 2.8 f-*. The results are consistent 
with a repulsive-core model of the deuteron. 
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DETERMINATION OF THE FORWARD SCAT- 
TERING AMPLITUDE FROM THE OPTICAL 
MODEL POTENTIAL IN THE HIGH-ENERGY 
LIMIT. R. Lipperheide and D. S. Saxon, De- 
partment of Physics, University of California, 
Los Angeles, California (Received June 6, 1960). 


A relation between the two-particle forward 
scattering amplitude, f,(0), and the volume in- 
tegral over the optical potential, Vp), is given 
in the high-energy limit, which is valid at least 
to second order in the two-body interactions, 
neglecting correlation and exchange terms. The 
relation is 


2nh? 
m 





Jv,@odz=-2 ay (0), 

where A is the number of nucleons in the target, 
and m is the reduced mass of the two-particle 
system. The connection between the shape of 
the optical model potential and the density dis- 
tribution of the target nucleus is discussed. 


SEARCH FOR S=-2 NEGATIVE HEAVY MESONS 
IN NUCLEAR EMULSIONS. Y. Eisenberg and 

M. Friedmann, The Weizmann Institute of Sci- 
ence, Rehovoth, Israel, and G. Alexander and 

D. Kessler, Israel Atomic Energy Commission 
Establishment, Rehovoth, Israel (Received 

June 24, 1960). 


A search for the D™ meson was conducted in an 
emulsion stack exposed to the 300-Mev/c K™ 
beam at Berkeley. The stars initiated by the 
beam particles were investigated for possible 
emission of two strangeness-carrying particles, 
and in addition mass measurements were done 
on the short-range tracks. In about 10000 events 
examined, no possible S=-2 particle has been 
observed. Hence the upper limit for the D~- 
meson occurrence in the present experimental 
conditions is ~1/100 from the double strangeness 
observation and ~ 1/6000 from direct mass deter- 
minations. 


MULTIPLE MESON PRODUCTION IN NUCLEON- 
ANTINUCLEON ANNIHILATIONS. LeRoy F. 
Cook, Jr.” and Joseph V. Lepore, Lawrence 
Radiation Laboratory, University of California, 
Berkeley, California (Received April 6, 1960). 


A two-parameter model is proposed for treat- 
ing complicated production problems in a rela- 
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tively simple way. It is assumed that the inter- 
action may be characterized by a range of inter- 
action and by a coupling strength. After the mode] 
is developed, it is applied to the problem of 
pion production in N-N annihilations. The two 
parameters are fixed by the experimental data 
for the multiplicity and energy spectra. It is 
found that all the data can be satisfied if one 
chooses the radius of interaction to be one pion 
Compton wavelength. Under certain restrictions 
the model reduces to the Fermi model. 





* 
Present address: Palmer Physical Laboratory, 
Princeton University, Princeton, New Jersey. 


INTERNAL PAIRS FOLLOWING 1” CAPTURE 
IN HYDROGEN. M. Derrick, J. G. Fetkovich, 
T. H. Fields, and J. Deahl, Carnegie Institute 
of Technology, Pittsburgh, Pennsylvania (Re- 
ceived June 17, 1960). 


Internal electron pairs from the two reactions 
(1) 77 +p—n+7°=n+y+et+e~ and (2) 1~+p—n 
+e+e~ have been studied in a hydrogen bubble 
chamber. 2184 cases were seen. A geometrical | 
cutoff selected 1523 of these as suitable for mo- 
mentum measurement. By an analysis of the 
momentum spectrum the Panofsky ratio was 
measured to be 1.51+0.10. The total intensity, 
momentum partition within the pairs, and dis- 
tribution in virtual photon mass are in essential 
agreement with the theoretical predictions of 
Kroll and Wada as recently extended by Joseph. 


RELATIVE IONIZATION OF PROTONS NEAR 
THE END OF THEIR RANGE. Hans Bichsel, 
Department of Physics, University of Southern 
California, Los Angeles, California (Received 
June 13, 1960). 





The influence of multiple scattering on strag- 
gling and consequently on the ionization produced 
by protons close to the end of their range is 
derived. The effect due to nuclear reactions is 
shown to be small at 340 Mev. The observed 
data at this energy are satisfactorily accounted 
for. 








PION-NUCLEON SCATTERING IN THE MAN- 

DELSTAM REPRESENTATION. S. C. Frautschi, 
Department of Physics, University of California, 
Berkeley, California, and J. D. Walecka, Insti- 
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tute of Theoretical Physics, Department of Phys- 
ics, Stanford University, Stanford, California 
(Received June 6, 1960). 


Using the analytic properties of partial wave 
scattering amplitudes, as derived from Mandel- 
stam’s representation, we have studied the 
J=3/2, P state of the pion-nucleon system. The 
method used is covariant, it incorporates uni- 
tarity, and the effect of a possible pion-pion re- 
sonance has been investigated. Using the “single 
nucleon term” and the low-energy scattering 
properties of the “crossed states,” we obtain a 
resonance in the J=3/2, T=3/2 pion-nucleon 
state without the aid of a cutoff. We have also 
investigated the scattering in the T=1/2 state. 
The pion-pion resonance appears to have only a 
very small effect in the T=3/2 state, whereas 
in the T=1/2 state it increases the phase shift 
by a factor of 2. 

The resonance obtained in the T=3/2 state oc- 
curs at too low an energy. There are several 
factors which may account for this: We have not 
been able to include fully the contributions from 
crossed states, and we have not systematically 
included inelastic scattering. 


GENERALIZED ANGULAR MOMENTUM IN 
MANY-BODY COLLISIONS. Felix T. Smith, 
Stanford Research Institute, Menlo Park, Cali- 
fornia (Received May 2, 1960). 


With short-range forces, initial and final states 
in a classical 3-body collision are straight-line 
trajectories into and out of a region where all 
three particles are close together at the same 
time. Using six coordinates, three describing 
the relative position of a pair of particles, and 
three the relative position of the third particle 
and the center of mass of the pair, the condition 
for simultaneous togetherness can be expressed 
with the help of the 66 grand angular momentum 





tensor, A, whose components are Ajj = (m;/mj;)” 
Xx5 Dj - (m;/m;)”*x ;p;. For a close 3-body col- 
lision, A* = 3)"; Ajj” must be small. A* com- 
mutes with the ordinary angular momentum oper- 
ators and with the kinetic energy; its eigenvalues 
are \(A+4) h?, with integral A, and its eigenfunc- 
tions hyperspherical harmonics. Initial and final 
3-body states can be described quantally by the 
total energy E, A?, and a commuting set of ordi- 
nary angular momenta; this description has the 
Same relation to a momentum representation as 
the ordinary angular momentum analysis has for 


a 2-body collision. A collision of (N+1) particles 
can be described by using a hierarchy of operators 
An* (2 <n <N); their eigenvalues are Ay(Ay, + 3 - 2) 
xh?, 


REMARKS ON SCHRODINGER’S MODEL OF DE 
SITTER SPACE. W. Rindler, Department of 
Mathematics, Cornell University, Ithaca, New 
York (Received June 6, 1960). 


Schrodinger in his book “Expanding Universes” 
developed a “reduced model” of de Sitter space- 
time consisting of a certain hyperboloid of one 
sheet. He showed, inter alia, that timelike sec- 
tions of the hyperboloid by planes through its 
center correspond to free paths in de Sitter space- 
time. The main purpose of the present note is 
to show that timelike sections by arbitrary planes 
correspond to paths of uniformly accelerated 
particles, and then to deduce some simple pro- 
perties of such paths. 


ABSENCE OF BOUND STATES IN A GRAVITA- 
TIONAL FIELD. Asher Peres, Department of 
Physics, Israel Institute of Technology, Haifa, 
Israel (Received June 30, 1960). 


The Klein-Gordon and Dirac equations for a 
particle in the gravitational field of a point mass 
are investigated. It is shown that the geometrical 
properties of the Schwarzschild metric prevent 
the normalization of any bound state solutions. 


QUENCHING OF MAGNETIC MOMENTS IN 
NUCLEI. S. D. Drell and J. D. Walecka, De- 
partment of Physics and Institute of Theoretical 
Physics, Stanford University, Stanford, Cali- 
fornia (Received June 20, 1960). 


Starting from the premise that with modern 
dispersion-theoretical techniques one has a 
reliable method for calculating the anomalous 
magnetic moment of a nucleon, we have calcu- 
lated the modification or “quenching” of this 
moment for a nucleon in nuclear matter. The 
effect we consider here is due to the fact that 
nucleons are not allowed by the exclusion prin- 
ciple to recoil into states already occupied by 
other nucleons in the nucleus. The actual tech- 
nique we have used in our calculation is to sum 
all the Feynman diagrams that are included in 
the dispersion theory calculation of the single- 
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